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Conventional near-shore large-scale hydrographic surveys 
use only two lines of position (LOPs) for position fixing. 
Previous works have proved that accuracies needed to meet 
the present International Hydrographic Organization 
standards are frequently not achieved using conventional 
surveying methods. The concept of using multiple lines of 
position (MLOP) adjusted by the least squares method was 
described. Actual field measurements acquired in the autumn 
of 1984 were processed to ascertain the increase in accuracy 
using MLOP versus conventional two LOPs on each hydrographic 
position. Recommendations are to use a four-range fully 
automated posit ion- fixing method to increase production, 
improve data quality, and have better control of the survey 
operations . 
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THESIS DISCLAIMER 



The reader is cautioned that computer programs developed 
in this research may not have been exercised for all cases 
of interest. While every effort has been made, within the 
time available, to ensure that the programs are free of 
computational and logic errors, they cannot be considered 
validated. Any application of these programs without 
additional verification is at the risk of the user. 



4 



DU 



\ 



NAVAL 

MONTEREY, Ca^i 



sonlVisi aoy^o-d002 



TABLE OF CONTENTS 



I. INTRODUCTION ...... 10 

A. BACKGROUND 10 

B. HYDROGRAPHIC STANDARDS OF ACCURACY AND 

METHODS ............ 12 

C. OBJECTIVES 14 

II. PROBLEM DEFINITION 16 

A. LARGE-SCALE SURVEY POSITIONING METHODS . .-. . 16 

1. Range-Range Positioning Method 17 

2. Azimuth- Azimuth Positioning Method .... 20 

3. Hybrid Positioning Methods 22 

B. ERROR THEORY AND ERROR PROPAGATION 24 

1. Types of Errors 24 

2. Accuracy of Large-Scale Survey Methods . . 29 

3. Error Propagation 34 

III. DATA ACQUISITION PROCEDURES AND SYSTEM 

PERFORMANCE 41 

A. DATA ACQUISITION PROCEDURES 41 

B. DATA ACQUISITION SYSTEM 46 

C. PERFORMANCE OF DATA ACQUISITION SYSTEM .... 46 



5 



IV. MATHEMATICAL MODEL 49 

A. LEAST SQUARES MODEL 49 

1. Unweighted Equations 49 

2. Weighted Equations 52 

B. PRECISION OF ADJUSTED VALUES 54 

V. DATA ANALYSIS AND RESULTS 57 

A. DATA PROCESSING PROCEDURES ..... 57 

B. DATA ANALYSIS 58 

1. Range Measurement Methods . 62 

2. Azimuth Measurement Methods 65 

3. Hybrid Methods • 69 

C. COST COMPARISONS OF POSITIONING METHODS ... 73 

D. OPTIMIZED SYSTEM 77 

VI. CONCLUSIONS 80 

A. MLOP FOR ON-LINE CALIBRATION OF RANGING 

SYSTEMS 80 

B. FUTURE ASPECTS AND CONSIDERATIONS 82 

APPENDIX A: FORTRAN PROGRAMS FOR DATA PROCESSING .... 84 

APPENDIX B: FORTRAN PROGRAM FOR LEAST SQUARES 

POSITION FROM MULTIPLE RANGES 87 

APPENDIX C: FOUR-RANGE MLOP FIX ACCURACIES 96 

APPENDIX D: FOUR- AZIMUTH MLOP FIX ACCURACIES 99 



6 



APPENDIX E: FORTRAN PROGRAM FOR LEAST SQUARES 

POSITION FROM HYBRID METHODS 103 

APPENDIX F: TWO-RANGE AND TWO-AZIMUTH MLOP FIX 

ACCURACIES 116 

LIST OF REFERENCES 119 

INITIAL DISTRIBUTION LIST 121 



7 



LIST OF TABLES 



I EQUIPMENT LIST 43 

II GEODETIC CONTROL POSITIONS 44 

III LOP COMBINATIONS FOR TWO HYBRID METHODS 60 

IV DATA FILE COMBINATIONS 61 

V POSITIONAL ACCURACIES FOR RANGE LOPS ...... 63 

VI POSITIONAL ACCURACIES FOR AZIMUTH LOPS 67 

VII POSITIONAL ACCURACIES FOR HYBRID LOPS 71 

VIII COMPONENT EQUIPMENT COSTS 74 

IX COST COMPARISON OF EQUIPMENT 75 



8 



LIST OF FIGURES 



2.1 Range-Range Positioning Method 18 

2.2 Azimuth-Azimuth Positioning Method 20 

2.3 Range-Azimuth Positioning Method 23 

2.4 Conversion of Angular Measure to Arc Length .... 33 

3 . 1 Survey Area 42 

3.2 Data Acquisition System 47 



9 



I. INTRODUCTION 



A . BACKGROUND 

A hydrographic survey is a compilation of numerous types 
of data which are acquired to ultimately produce a nautical 
chart. The two most important aspects of a hydrographic 
survey are (1) the depth sounding of the seabed and (2) the 
simultaneous acquisition of a horizontal position which is 
attached to the sounding. If either of these two aspects 
are in error then the overall accuracy of the survey is 
lost. This thesis will address the latter aspect. Ingham 
[1984, p. 14] states 

A major element of the hydrographic survey is the 
dynamic positional control of the survey vessel. 

Indeed, the great majority of the surveyor's tasks 
offshore today require positioning services alone. 

The accurate determination of a position at sea has 
always been a painstaking effort for the hydrographic 
surveyor. First, a survey of the coast adjoining the 
offshore survey area must be made to establish an accurate 
position of each shore control site. This shore control 
must then be "extended" offshore by a variety of means and 
systems depending on the scale and purpose of the survey. 

The position of a survey vessel is either determined by a 
visual method, such as a three-point sextant fix to visual 
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signals onshore, or by a radio navigation method. In the 
very near future, the Global Positioning System (GPS) 
promises to completely revise the method of offshore 
navigation and surveying. Consistent differential 
positional accuracies of about 10 m, or less, are expected 
in real time once the system becomes fully operational thus 
allowing the hydrographer the freedom of surveying as far 
offshore as the situation demands. 

This thesis will address various ways to increase the 
positional accuracy of nearshore surveys. Nearshore areas 
(i.e., channels, harbors, and bays) pose the most hazard to 
the mariner and are the areas where accurate position 
determination of the bottom features is essential. 

Nearshore surveys of 1:5,000 scale, and larger, require 
intensive and time consuming delineation of shoreline 
features and bottom topography with positional accuracies of 
5 m, or less. Presently, GPS will not provide these 
accuracies in a real-time mode. Conventional hydrographic 
survey methods are the only way to determine the position of 
a survey vessel in real time. Mobley stated [Saxena, 1980, 
p. A-4] , 

Positioning systems for precise navigation for survey 
operations are barely adequate for the scales of surveys 
being accomplished today. A lot of time and effort is 
lost due to equipment down time, moving equipment and 
the use of only 2 lines of position . . . which the 

present National Ocean Survey (NOS) DAS can only accept. 
New and future replacement systems providing redundant 
LOPs are a must if higher accuracies are desired 
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offshore. Global Navigation may be the answer but will 
be 5 to 10 years before readily available to all users. 

If conventional systems and methods are to be used, a 
redundancy in the measurements is necessary to obtain a 
higher order of accuracy than the present system of 
acquiring two lines of position (LOPs) on each fix. This 
thesis will investigate the application of multiple lines of 
position (MLOP), adjusted by the least squares method, for 
hydrographic position control. 

B. HYDROGRAPHIC STANDARDS OF ACCURACY AND METHODS 
The National Ocean Service (NOS) is the agency 
responsible for nautical charting in the United States and 
has documented many standards and procedures for the 
acquisition of data during a hydrographic survey. With 
respect to the position of the survey vessel, the 
Hydrographic Manual [Umbach, 1976, p. 4-14]‘ states 

Hydrographic position fixes on sounding lines are almost 
always determined by the intersection of two lines of 
position; dead- reckoning positions based on course and 
speed provide an additional internal check. 

Also Ingham [1984, p. 15] states in his chapter on 
Posit ion- fixing Afloat 



For any fix one requires at least two lines of position 
. . . obtained from readily identifiable sources. These 

should intersect at the largest angle possible for good 
fix-geometry, and have the smallest possible standard 
deviation, resulting in the smallest area of 
uncertainty . 
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NOS requires that for surveys of 1:5,000 scale, and 
larger, at least one of the LOPs must be from a visual 
method. Thus, a standard way of acquiring positional fixes 
on large-scale surveys has been the range-azimuth 
(rho-theta) method as addressed by Waltz [1983] . One of the 
reasons this hybrid method is popular is because the two 
LOPs will always intersect at 90°, the best intersection 
possible for good accuracy. 

Unlike mariners, the hydrographer has access to more 
sophisticated and more accurate positioning systems. But as 
Riemersma [1980, p. 10] observes 

There is no mariner or surveyor in the world who in 
conventional navigation does not use a check bearing 
just to make sure. Why then does one, with 
radiopositioning [sic] systems, use only two LOPs? 

By using only two LOPs there exists a level of uncertainty 
with each position which can only be resolved during the 
post processing phase of the survey [Perugini, 1984]. In 
the past, two LOPs were the most feasible and sometimes the 
only way to determine the position of the survey vessel. 

With the rapid advance in electronic technology, the 
acquisition of multiple-range LOPs are as common today as 
two- range LOPs were 20 years ago. However, due to cost 
limitations, or possibly the aversion of changing tried and 
true operational procedures , hydrographic surveys are still 
usually controlled with only two LOPs. The acquisition of 
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only two LOPs on each fix will change in the very near 
future for NOS as they are in the process of implementing a 
specification from the 1982 International Hydrographic 
Organization (IHO) [Anonymous, 1982, p. 5] stating 

It is desireable that whenever positions are determined 
by the intersection of lines of position, three such 
lines be used. The angle between any pair should not be 
less than 30° . 

With MLOP the hydrographer will have redundant 
observations on each fix. These redundant measurements may 
be adjusted, as in normal land surveying, by the least 
squares method to acquire the "best" fit. Applying the 
least squares method to hydrographic position adjustment 
allows the hydrographer to statistically determine errors in 
real time if computer programs are operated aboard the 
survey vessel [Silva, 1982, p. 100], 

C. OBJECTIVES 

Will the use of MLOP increase the accuracy of the fix 
over the present method of two LOPs? This question has been 
addressed by both Kaplan [1980] and Silva [1982] who found 
that MLOP do improve the accuracy of the survey. However, 
the method was applied to theoretical data only. Therefore, 
one objective of this thesis was to acquire data for a 
portion of an actual hydrographic survey with a minimum of 
four LOPs on each position fix. A model was constructed to 
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determine the least squares position for each fix from the 
MLOP. The precision of each MLOP fix was compared to the 
precision of a fix obtained by two LOPs to determine if a 
major difference did exist between the two methods. 

To obtain the additional LOPs, a significant increase in 
the cost of operations can be expected and must be 
justified. More LOPs will require more horizontal control 
for shore reference points. Additional survey equipment, 
electronic range measurement units or azimuth measurement 
units, or both, will be necessary to acquire the four LOPs. 
Initial equipment costs could possibly double at the onset 
of this method of surveying. Consequently, all of the costs 
associated with the MLOP method must be weighted against the 
optimal position accuracy requirement to determine if the 
additional accuracy is feasible within the scope of a normal 
hydrographic survey. 

A final objective of this thesis was to determine an 
optimal system which is the most cost effective as well as 
one which will meet the new IHO standards for position 
fixing of soundings . 
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II. PROBLEM DEFINITION 



The geometry of a standard position fix using two LOPs 
must be understood before the various MLOP methods may be 
evaluated. Therefore, the configuration of large-scale 
positioning methods, the errors unique to the equipment 
used, and the propagation of these errors into the offshore 
position warrant discussion in this chapter. 

A. LARGE-SCALE SURVEY POSITIONING METHODS 

The intersection of two lines defines a point in space. - 
In navigation and offshore surveying each of the two lines, 
or LOPs, can be generated via an electronic or visual 
method. Traditional hydrography requires that the two LOPs 
generated intersect at an angle greater than 30° and less 
than 150°. 

Notations used for computations are: 

a) X coordinate positive to the east, Y coordinate 
positive to the north 

b) Occupied Shore Stations; numerals 1, 2, . . . 

c) Unoccupied Shore Station used for Initial Azimuths; 
numeral 0 

d) Inverse Distance between shore stations; D 

e) Ranges to vessel, point P; R-^, R£ , . . . 

f) Computed horizontal angle; a 

g) Measured horizontal angles; a^» ct2 > • • • 

h) Azimuths of lines; 



a l , 2 ’ a l,P> ‘ * 
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1 . Range-Range Positioning Method 



An electronic range-range positioning method on 
large-scale surveys is generally accomplished with 
short-range (line-of-sight ) equipment. At these short 
distances the ranges to the vessel are determined by 
measuring the time of transit or difference in transit times 
of an electromagnetic pulse [Davis et al . , 1981]. The 
transit-time measurement starts with an interrogation pulse 
from a master transponder onboard the vessel which is 
received by a remote transponder erected over a known 
control point and returned to the master transponder. The 
total transit time is then converted to a distance by 
assuming an average value for the propagation velocity of 
the signal. The LOPs generated from a system such as Del 
Norte Trisponder (used for this thesis) are concentric 
circles from the shore stations. The Del Norte system 
operates at a frequency of 9 GHz with a pulse repetition 
interval (PRI) which can be selected from 304 to 998 
microseconds. Thirty-two consecutive valid time-difference 
measurements are processed to provide each range output 
[Ingham, 1984, pp . 37-38]. 

The geometry of a range-range position fix using two 
shore stations is graphically illustrated in Figure 2.1. 
Basic trigonometric relationships are employed to determine 
the position of the survey vessel at point P using the 
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Figure 2.1 Range-Range Positioning Method. 

measured ranges R-^ and R£ from shore control stations 1 and 
2, respectively. 

The computation of the coordinates of the survey 
vessel, Xp and Yp, at point P is as follows [Davis et al., 
1981, p. 876] : 

D = J (X 1 -X 2 ) 2 + (YrY 2 ) 2 (2-1) 

where X-p Y-^ are the coordinates of shore station 1 and X£ , 
Y£ are the coordinates of shore station 2. Since all three 
sides of the triangle are known, angle a can be computed 
from the law of cosines equation 
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a = arc cos [(D 2 * Rj 2 - R 2 2 )/(2 R x D)] 



( 2 . 2 ) 



The azimuth of the line from station 1 to station 2 can be 
determined by the inverse computation equation 



a l,2 = arc tan C( x 2~ x l) / ( y 2" y 1^ 

The azimuth from shore station 1 to the vessel P is 


(2.3) 


a l,P = a l , 2 + a 


(2.4) 



If Xp and Yp are the known coordinates of the left-hand 
station when looking from the vessel toward the base line, 
then the coordinates of the vessel at P can be computed by 
the equations 



Xp = Rp sin ( cci 5 p ) + X 1 


(2.5) 



Yp = Rp cos ( ap^p ) + Yp 


(2.6) 


and the angle of intersection of the two range LOPs 
computed by the equation 


can be 


J3 = 180° - arc cos[(D 2 - R-^ - R 2 2 ) / (2 Rp R 2 )] 


(2.7) 



Determining the position of a vessel using multiple 
ranges from three, or more, shore stations will involve 
additional modeling since a redundancy of LOPs are 
generated. Chapter IV will address this problem in detail. 
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2 . Azimuth-Azimuth Positioning Method 



An azimuth- azimuth positioning method is a direct 
extension of the standard land surveying practice of azimuth 
intersection. Ingham [1984, p. 16] states 

This [method] is usually not to be recommended, since 
complex arrangements are required to ensure that both 
angles are measured at the instant of the fix marked on 
the echo-sounder record, and to communicate the fix data 
to the boat. 




Figure 2.2 Azimuth- Azimuth Positioning Method. 
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Although true in some instances, good radio communications 
will alleviate the problem. Another option is to use a 
flagman aboard the survey vessel to coordinate the 
activities . 



theodolites, set up over known control stations, which are 
sighted upon the survey vessel as it moves along the 
sounding line. When done correctly, the method produces a 
very accurate position at the intersection of the visual 
LOPs. The azimuth to the vessel is obtained by initializing 
on a known control station and then measuring the clockwise 
angle to the survey vessel (Figure 2.2). If theodolites 
occupy the known geodetic shore stations 1 and 2, and 
station 0 is the station which is initialed upon, then the 
coordinates of point P, Xp, and Yp can be determined by the 
point-slope equations [Davis et al . , 1981, p. 371]. 



The azimuth-azimuth positioning method uses two 



Yp - Y-^ + (Xp - X x ) cot ( a^p ) 



( 2 . 8 ) 



and 



Y p = Y 2 + (X p - X 2 ) co t ( a 2 , P ) 



(2.9) 



Subtraction of Equation 2.8 from Equation 2.9 yields 




[(Y X - Y 2 ) - X-l cot ( a 1>P ) 

+ X 2 co t ( a 2 s p )] / 

[cot ( a 2jP ) - co t ( ax jP )] 



( 2 . 10 ) 
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and 



Y P = Y 2 + (X p - X 2 ) cot( a 2)P ) (2.11) 

Figure 2.2 shows 

<*1,P = a l,0 + a l (2.12) 

and 

a 2,P = a 2 , 0 + a 2 (2.13) 

The use of MLOP obtained from multiple azimuths to 
determine the position of the vessel is discussed in detail 
by both Kaplan [1980] and Silva [1982] . 

3 . Hybrid Positioning Methods 

A hybrid positioning method (i.e., range-azimuth 
method) is a frequently used positioning method on 
near-shore large-scale surveys. The method provides an 
excellent geometric determination of the vessel's position 
since the two LOPs always intersect at 90°. One LOP will be 
generated from a measured electronic range to the vessel. 

The other LOP is the observed azimuth to the survey vessel. 

A concentric configuration of the ranging and angle 
measuring devices is commonly used on large-scale surveys. 

A concentric configuration is not always possible due to 
terrain and the location of the control stations. In that 
case, one of the measurements must be made from an eccentric 
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p 




Figure 2.3 Range-Azimuth Positioning Method. 

station IE (Figure 2.3). In practice station 1 and station 
IE are within about 2 m of each other and are occupied by a 
theodolite and a remote ranging unit, respectively. 

Knowing the coordinates of stations 1 and IE, the 
coordinates of point P may be computed by the equations 
[Wallace, 1971, pp . 46-48] 

Xp = + R]_ sin( Ct^E^p ) (2.14) 

Yp = Y^£ + R^ cos ( ct]_E 5 p ) (2.15) 
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where p is derived from the equation 



sin( ~ a^p ) - [1 / R]J (2.16) 

cos ( a^p ) _ ( ) sin ( a^ p )] 

which reduces to 

a^E ^ p = ctp s p + arc sin ([(X^-X^) cos( a^p ) (2.17) 

- (Y x -Y 1e ) sin( a 1>P )] / R x ) 

B . ERROR THEORY AND ERROR PROPAGATION 

The present NOS Hydrographic Manual [Umbach, 1976, p. 
1-4] states that positional error on any fix shall seldom 
exceed 1.5 mm at the scale of the survey. Normally 1.0 mm 
is allocated for navigational system error and 0.5 mm for 
plotting error. Thus, for a 1 : 5 , 000- scale survey the total 
allowable navigational error is 5 m. The rated accuracy of 
most positioning systems used approaches this allowable 
error. Errors must be identified and dealt with 
accordingly. 

1 . Types of Errors 

Types of errors as applied to hydrographic surveying 
are described [Davis et al., 1981; Greenwalt and Schultz, 
1962; Heinzen, 1977; Kaplan, 1980; Mikhail, 1976; Perugini, 
1984; Waltz, 1983]. 

a. Blunders 

Blunders, also known as gross errors or 
mistakes, are mainly due to equipment malfunctioning or 
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observer carelessness. Blunders usually are of large 
magnitude in comparison to other errors . Unlike land 
surveying, hydrographic surveying is a dynamic situation and 
the existence of blunders in field data is commonplace. 

When using a theodolite for position control of the vessel, 
the observer rarely has time to check the pointing to the 
moving target (survey vessel). Also the acquisition of 
multiple readings and repeating measurements independently 
to check for consistency, which are standard procedures in 
land surveys, are not feasible. On any particular fix, an 
observer has just one chance of an accurate sight to the 
vessel. If any doubt exists as to the quality of the 
sighting, the observation should be rejected. 

Blunders can occur when using electronic 
positioning systems due to malfunctioning of the equipment. 

A problem inherent with short-range systems, such as Del 
Norte Trisponder, is multipath effects (reflected wave from 
the ocean surface as well as the direct transmission) 
[Munson, 1977, p. 4]. Fading, dropout of the ranging 
signal, and reflected signals are commonly observed during 
field work due to the omnidirectional antenna aboard the 
survey vessel. 

Detection of blunders is one of the daily 
functions of a hydrographer . Automated data acquisition 
systems have the capability of "flagging” blunders on-line 
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through software engineering. During nonautomated 
surveying, a plot of the position fixes is maintained to 
check the accuracy of each fix as data are acquired. Most 
blunders occur during the nonautomated phase of data 
acquisition due to the increased amount of hand logging. 

Blunders can be minimized by using automatic 
electronic systems for data acquisition, such as the NOS 
HYDROPLOT Data Acquisition System (DAS) or the Racal Decca 
AUTOCARTA II system. Also, all blunders must be detected 
and resolved prior to the final processing of hydrographic 
data. Normally, repeated checking of the data will minimize 
the amount of blunders. 

b. Systematic Errors 

Systematic errors follow a definite pattern and 
are generally constant in magnitude and sign throughout a 
series of observations. The system causing the pattern may 
be dependent on the instrument, or atmospheric effects 
[Mikhail, 1976, p. 67]. In hydrographic surveying, 
systematic errors must be determined and data corrected 
prior to the final processing adjustment of all the survey 
records . 

Determination of systematic errors in 
hydrography is accomplished through the calibration of the 
surveying instruments. Calibration is the process of 
comparing the observed instrument value with a known 
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standard. The difference between the observed and known 
values is the total systematic error present in the system. 
Once the systematic errors are determined, a "corrector" is 
applied to the observed values which yields the "true" value 
of the observations. A corrector is equal in magnitude but 
opposite in sign to the systematic error. 

A systematic error existed within the Del Norte 
Trisponder microwave ranging equipment used as the 
positioning system for data acquired for this thesis. The 
standard field procedure for this equipment was used to 
compare ranges obtained from the Trisponders to a measured 
geodetic base line. Differences were recorded and applied 
during the post processing phase of the survey as correctors 
to the ranges. Daily calibration checks were made while the 
vessel was underway to ensure that no major changes had 
developed, since noticable calibration drifts can occur with 
time and variation in range [Munson, 1977, p. 4]. Base-line 
calibrations of short-range systems are normally performed 
before, during, and after the completion of data acquisition 
on a hydrographic survey. 

Numerous systematic errors arise when using 
theodolites for range-azimuth and azimuth- azimuth 
positioning on large-scale surveys. Proper adjustment and 
alignment of the instrument will eliminate a number of 
errors. Mikhail [1976, p. 71] lists some of the common 



27 



sources of systematic errors associated with the use of 
theodolites. Initial pointing and eccentricity are two 
additional sources of errors introduced in hydrographic 
surveys . 

An initial pointing must be made to a known 
geodetic station since the observer is actually measuring a 
direction to the survey vessel. Normal practice in the 
field is to adjust the plate setting to a number other than 
000 o 00'00", usually about 000°00'10 M . This value must be 
recorded in the record volume of the survey and applied as a 
corrector during processing. Should this initial corrector 
be ignored, or applied in the wrong fashion, a systematic 
error will be present in all of the observations. 

Eccentricity frequently occurs when using a 
range-azimuth positioning method as it is very difficult to 
erect a remote ranging unit (i.e., a Trisponder) and a 
theodolite directly over a horizontal control station. 

Unless two different height tripods are erected over the 
station, or a platform built to hold both instruments, an 
eccentric horizontal station must be established to locate 
one of the survey instruments. If this offset position is 
not accounted for, then a systematic error (varying in 
magnitude with the range and direction to the vessel and the 
eccentricity of the station) will result. During the data 
acquisition phase for this thesis, the ranging unit used for 
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observing the range to the vessel was erected over an 
eccentric station at three of the four horizontal control 
stations . 

c. Random Errors 

Random errors account for the remaining 
variation in measurements after all blunders have been 
discovered and removed, and the observations corrected for 
the known systematic errors . Random errors cannot be 
modeled mathematically but must be modeled under the laws of 
probability. The random error is a random variable which 
can take on a number of possible values depending on the 
probability involved. 

They result from accidental and unknown combinations of 
causes beyond the control of the observer. Random 
errors are characterized by: (1) variation in sign -- 

positive and negative errors occuring with equal 
frequency, (2) small errors occurring more frequently 
than large errors, and (3) extremely large errors rarely 
occuring [Greenwalt and Schultz, 1962,. p. 2]. 

Statistical treatment of random errors relating 
to hydrographic applications has been given by Waltz [1983] 
and Perugini [1984] . 

2 . Accuracy of Large-Scale Survey Methods 

Two LOPs from completely automated electronic 
systems can rarely meet the required accuracy for a 
1: 5 ,000-scale survey due to the standard error associated 
with each LOP [Perugini, 1984]. Common survey practice has 
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been to use hybrid methods, such as range-azimuth, to reduce 
this error by introducing an accurate direction measured 
with a theodolite which will have a 90° intersection with 
the range LOP. Also completely visual methods can be 
employed, such as azimuth- azimuth or three-point sextant 
fixes, for the position determination of the survey vessel. 
All of the methods employed contain unique standard errors 
associated with the equipment utilized. This section 
discusses the rated accuracy of the equipment used for data 
acquisition for this thesis. 

a. Errors in Ranges 

The total positional accuracy "shall seldom 
exceed 1.5 mm at the scale of the survey" [Umbach, 1976, p. 
1-4]. Munson [1977, p. 2] has equated "seldom" to be 90% of 
the time which results in a 1.645-sigma value of 4.5 m on a 
1 : 5 , 000- scale survey. He further states 

Positioning system accuracies are most commonly stated 
in terms of error along a line of position. S ince the 
total position error will always be at least ^ j 2 
greater than this, an acceptable positioning system must 
measure a lin e of position to at least a la level of 
4.5 m / ^J ~~ 2 = 3.2 m for a 1:5,000 survey. 

The Del Norte Trisponder system used for this 
thesis has a manufacturer stated accuracy of ±1 m. Recent 
stability testing by NOS on the Model 520 Digital Distance 
Measuring Unit (DDMU) has shown errors in the order of 3 m 
over various ranges [Whitsell and Berstis, 1983, p. 4]. The 
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study also found that the range error was a function of low 
signal strength which could be simulated by signal 
attenuation. This is a major problem because 



There is no indication in the Model 520 DDMU display of 
data output when low signal strength conditions are 
encountered. Reduced signal attenuation will occur 
during hydrographic survey operations due to dynamic 
movement of the master antenna on the survey launch, 
inclement weather conditions, and the presence of 
multipath propagation zones. [Whitsell and Berstis, 
1983, p. 5] 



A 3-m standard error will be used for all range 
LOPs acquired during the field work for this thesis with a 
Model 540 DDMU. This unit is similar to the Model 520 DDMU, 
the only difference being that . the 540 has a four-range 
digital display versus a two-range digital display on the 
520. A fixed value over all ranges was chosen rather than 
one which varies with distance, as is normally defined for 
precise Electronic Distance Measuring (EDM) equipment in 
land surveying, to simulate common field practice. As a 
result, the error of each LOP at a short range (1,000 m) 
will have a higher bias than a LOP at a longer range (10,000 
m) due to the fixed error. 

b. Errors in Azimuths 

One of the objectives of the work done by Waltz 
[1983] was to determine the pointing error of the Wild T-2 
theodolite when used for range- azimuth hydrographic control. 
His rigorous treatment of the subject found that an 
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estimated linear error for the pointing of a T-2 was 1.3 m. 
Actually the pointing error is a function of the angular 
resolution of the instrument plus the ability of the 
observer to track the sounding vessel. Both are a function 
of the distance to the vessel. For azimuth- azimuth 
positioning methods, the errors in positioning actually 
depend upon the distance, the angular resolution, and the 
angle of intersection of the LOPs [Heinzen, 1977, p. 55]. 

An important conclusion made by Waltz [1983, p. 

81] was 

There exists, for any angle measured with a T-2, a time 
lag of about one second between angle observations and 
any measurement made aboard the vessel, including both 
automatic and manually recorded depth and range data. 

There is then an associated position error for these 
measurements, the magnitude of which depends upon vessel 
speed, which was about two meters for the four knot 
speed used in this experiment. 

With this fact in mind, a standard error 
associated with azimuth LOPs was computed for the data 
acquired for this thesis. Hand plots of the position fixes 
show that the speed of the survey vessel was approximately 5 
knots (2.6 m/ s ) . Thus, the 1- second time lag would amount 
to as much as an additional 2.6 m of error in position due 
to the vessel velocity relative to the shore station. 
Consider an example which shows positional error as a 
function of the range and the angular resolution of the 
instrument (Figure 2.4). Using the standard conversion 
formula 
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Figure 2.4 Conversion of Angular Measure to Arc Length. 

a / p = b / R (2.18) 

where b is the arc swept by the vessel (2.6 m) , R is the 
range .to the vessel in meters, a is the angular error in 
seconds, and p is a conversion factor (206265"/rad) . if R 
is 1,000 m, then a = (206265” ) (2 . 6 m / 1,000 m) = 536" = 

8.9'. Likewise, if R is 2,000 m, then a = 4.5'; if R is 
3,000 m, then a = 3.0'. 

A total standard error for a LOP measured by an 
azimuth for this thesis was assumed to be 4 m. This value 
takes into account a 1.3-m pointing error plus the 2.6-m 
error due to the 1-second time lag at a vessel speed of 
approximately 5 knots. This was done to duplicate standard 
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processing techniques for azimuth data. By assuming a fixed 
error over all ranges a reverse situation holds for LOPs 
from azimuth measurements when compared to the similar 
situation for range measurements. At longer distances the 
measured azimuth will be more accurate than at shorter 
distances because the aspect of the vessel does not change 
as rapidly. By using a fixed error in position, the 
azimuths obtained at close range will have an underestimated 
error whereas those observed at longer ranges will have an 
error which is overestimated. 

3 . Error Propagation 

Since random errors will exist in the field 
measurements, then the errors in the computed positions of 
the survey vessel can be obtained by the law of propagation 
of errors. Error propagation is the evaluation of errors in 
the computed positions as functions of the errors in the 
field measurements. As stated previously, random errors 
follow statistical laws and can be dealt with through 
probability distributions of these errors. 

Let the function Y be defined by 



Y = f (X 1 ,X 2 ,X 3 , . . . ,X n ) (2.19) 

in which X-p X£ , X 3 , . . . , X n are uncorrelated 

measurements. Thus, the law of propagation of errors 
[Mikhail and Gracie, 1981, pp . 148-178] can be shown as 
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( 2 . 20 ) 



a Y 2 = ( 3Y / dX-L ) 2 (j x 2 ♦ ( 3Y / 5X 2 ) 2 c x 2 

♦ . . . ♦ ( 5Y / 3X n ) 2 ct x 2 

n 

Consider the problem of defining the position of a 
survey vessel at sea. The goal of the positioning methods 
previously discussed is to accurately determine the computed 
coordinates (Xp Yp) of the survey vessel. Each of the two 
coordinates are independently a function of the measurements 
used to compute the coordinates. For example, when using a 
range-range positioning system the coordinates are a 
function of the measured ranges (Rp and R 2 ) and the known 
station coordinates on shore such that, 

X p = f 1 (R 1 ,R 2 ,X 1 ,X 2 ) (2.21) 

Y p = f 2 (R 1 ,R 2 ,Y 1 ,Y 2 ) . (2.22) 

where Xp, Yp and X 2 , Y 2 are the known coordinates of the 
shore stations 1 and 2, respectively. The law of 
propagation of errors takes the following form 

a x 2 = ( axp/aR-L ) 2 a Rl 2 * ( ax p /aR 2 ) 2 CTR 2 (2.23) 

* ( axp/axp ) 2 ct Xi 2 + ( ax p /ax 2 ) 2 CTx 2 

a Y 2 = ( aYp/dRx ) 2 cr Rl 2 * ( 3Y p /5R 2 ) 2 CTr 2 (2.24) 

* ( aY p /3Y 1 ) 2 a Yi 2 * ( aY p /aY 2 ) 2 CTY 2 
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The computed position (Xp, Yp) from a range-range 
positioning system was given by Equations 2.1 through 2.7. 
Partial differentiation of the terms in Equations 2.23 and 
2.24 yields Equations 2.25 through 2.28 where 
u = (D 2 + R x 2 - R 2 2 ) / (2 R x D ) 
then 



ax p / 5R X = sin( a 1>P ) - [cos( a ljP )] (2.25) 

( 1 - u 2 J’ 1 / 2 [ (R p 2 * R 2 2 - D 2 ) / (2 R x )] 

ax p / 3R 2 = [ co s ( a ljP )] ( 1 - u 2 )' 1 I 1 (2.26) 

(R 2 / D) 

5Y p / 5R]_ = cos( ai,p ) + [sin( ap iP )] (2.27) 

( 1 - u 2 r 1 / 2 [(Rp 2 + R 2 2 - D 2 ) / (2 R p )] 

aY p / <5R 2 = [sin( a 1)P )] ( 1 - u 2 )‘ 1/2 (2.28) 

(- R 2 / D) 



Since the X and Y coordinates of the shore stations are 
constant then 



ax p /ax 1 = dX p /dX 2 = aYp/^Yp = 5Y p /aY 2 =0 (2.29) 

Previously, it was shown that the standard error of 
the Del Norte Trisponder system, <jp, has been determined to 
be 3 m, regardless of range, so 



aR l ” aR 2 = a R ~ 3 m 



(2.30) 
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Then, the standard error of each of the computed coordinates 
will be 



a X = a R 




(3x p / aRj^) 2 ♦ (ax p / aR 2 ) 2 



( 2 . 31 ) 



and 



(J Y = Or 




(<3y p / aRx) 2 + (aY p / ^r 2 ) 2 



( 2 . 32 ) 



For range- azimuth positioning systems, the 
coordinates of the vessel are a function of the measured 
range and azimuth on each fix plus the X, Y coordinates of 
the range and azimuth station 



X 



P = ^l( R l»ai,p> x lE» x l) 



( 2 . 33 ) 



and 



Y p = f 2 ( R l’ai,P’ Y lE’ Y l) 



( 2 . 34 ) 



so that the equation of the propagation of errors takes the 
form 



<*X 



2 _ 



a Y 



2 . 



( aXp/oRi ) 2 ct R i 2 * ( ax p /aa liP ) 2 a ai p 2 * (2.35) 

( ax p /3X 1E ) 2 a XlE 2 * ( 3Xp/3 Xl ) 2 a Xl 2 

( aYp/SRx ) 2 ct r 2 * ( aYp/3a liP ) 2 a 2 * (2.36) 

( 3Y p /5Y 1E ) 2 «j Y ie 2 ♦ ( 3Yp/3 Yl ) 2 a Yl 2 



37 



The position coordinates (Xp Yp) were given in Equations 
2.14 through 2.17. The value of p (Equation 2.17) is 

repeated here (Equation 2.37) for continuity. 

®1E,P = arc s i n ([(X^-X^jr) cos ( ai jP ) (2.37) 

~ (Yx"Y^g) sin( ax , P )] / Rx^ + a l,P 

Partial differentiation of the terms in Equations 2.35 and 
2.36 where 

u = (1/R l ) [(Xx - X 1E ) cos (ax ^ P ) - (Y L - Y 1E ) sin(ax jP )] 
yields 



5Xp / d^x ~ s in (aiE ? p ) ~ u [cos (ax E ? P ) ] 



(2.38) 



( 1 - u* ) 



2 a-1/2 



£X p / dax, p = [cos(a 1E)P )] {^x " ( 1 " u2 )' 1/2 (2.39) 



(X x ~ Xx E ) [sin(ax 5 p)] + (^x ~ ^1E^ cos (ax )P )} 



5Y p / 5^x _ cos (ax E ^ p ) + a [ s in ( a x E ^ P ) ] 



2 'k - 1/2 



(2.40) 



( 1 - u* ) 



5Yp / 5ax 5 p ~ [sin(ax E , p)l {[( 1 - u ) 



2 x-1/2 _ 



R x ] (2.41) 



( Y 1 ~ y 1 E ) C cos (a x 5 p) ] + (X x - X 1E ) sin(ax jP )} 
Since the X and Y coordinates of the range and azimuth shore 
stations are constants then 



ax p /ax 1E = axp/axp = 5Y p /dY 1E = dY ? /dY l = o 



(2.42) 
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Substitution into the propagation of error equations 
(Equations 2.35 and 2.36) results in 

<TX^ = ( ^Xp/dRp ^ + ^ ^ x p/ dap^p o^p^p^ (2.43) 

a Y 2 - ( aYp/^Rx ) 2 a R 2 + ( 3Y P /3ai ; p ) 2 a al>P 2 (2.44) 

and the standard error ax a nd cry of the X and Y coordinates 

is obtained from Equations 2.43 and 2.44. 

The standard error of position, ap > was used as a 
reference to classify the accuracy of the data acquired for 
this thesis. Rather than classify the accuracy using the 
standard error ellipse or confidence circles , the Qp value 
was chosen due to its ease of computation. The value of the 
standard position error [Saxena, 1972, p. 15] is defined as 
the square root of the sum of the squares of the propagation 
variances of the system. Mathematically 

CT p = J CT X 2 + °Y 2 ~ (2.45) 

p p 

where ax and cry are computed as shown previously for 
range-range LOPs and range- azimuth LOPs . The standard error 
of position, a s in cm, for azimuth LOPs [Mueller, 1979, p. 
62] is 

a s = ( 0.485 / sin p ) a a J Dp 2 + °2 2 (2.46) 



39 



where 0 is the angle of intersection of the two LOPs ; is 
the distance from point 1 to point P in km; T>2 is the 
distance from point 2 to point P in km; and <j is the 

vA 

angular error in seconds of arc. 

For the overdetermined case of MLOP, the <jp value 
(as derived by Equations 4.28, 4.29, and 4.30) will be 

cr x = a o V Q xx (2.47) 

ct y = a 0 yj Qyy (2.48) 

and 

tfp = V (?x 2 + ct Y 2 = a o V Qxx + Qyy (2.49) 

where <j 0 is the standard error of an observation of unit 
weight , Qxx and Qyy are the diagonal elements of the 
variance- covariance matrix. 

Thus a tag can be placed on every position fix 
acquired from whatever system, or combination of systems, 
used to determine the coordinates of the survey vessel. 
Using this standard position accuracy, a comparison between 
fixes determined by different systems and methods can be 
compared so that the most precise combination of ranges, 
azimuths, or hybrid methods can be analyzed. 
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III. DATA ACQUISITION PROCEDURES AND SYSTEM PERFORMANCE 



The field work required for this thesis involved 
conducting a hydrographic survey while using MLOP for 
sounding line position fixing. Two types of positioning 
data, ranges and azimuths, were acquired simultaneously from 
four shore control sites. The data were then analyzed using 
various combinations of ranges and azimuths to determine a 
more accurate positioning system. 

A. DATA ACQUISITION PROCEDURES 

Students in the NPS Hydrographic Sciences curriculum 
conducted a basic hydrographic survey of southern Monterey. 
Bay as part of a field experience requirement.. The survey 
was accomplished using procedures similar to those used by 
NOS for nearshore surveying. A chartered 36-foot, 
twin-engine Uniflite boat was the platform used for the 
survey. During the week of 26 November 1984 to 30 November- 
1984, the boat was positioned using four ranges and four 
azimuths from known third-order horizontal geodetic control 
stations on shore for each position fix of the vessel. The 
ranges to the vessel were determined from four Del Norte 
Trisponders and were automatically recorded by a Racal-Decca 
AUTOCARTA II data acquisition system. The four azimuths 
were observed with Wild T-2 theodolites. Data acquisition 
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Figure 3.1 Survey Area. 



equipment was supplied by Racal-Decca, NOAA, and NPS 
(Table I) . 



TABLE I 

EQUIPMENT LIST 



Aboard M/V Silver Prince 

Model 540 DDMU S/N 211 (Racal) 
Master Trisponder S/N 3014 (NOAA) 



Onshore Control Stations 



STATION 


REMOTE TRISPONDER 


T-2 


(S/N) 


SQUARE 1984 


Code 72 

S/N 2819 (Racal) 


14405 


(NPS) 


CONK 1984 


Code 74 

S/N 2822 (Racal) 


51642 


(NPS) 


USE MON 1978 


Code 7 6 

S/N 3004 (NOAA) 


30504 


(NPS) 


GEOCEIVER 1982 


Code 78 

S/N 2986 (NOAA) 


14482 


(NPS) 



Geodetic control for the survey consisted of third-order 
monumented stations and established eccentrics (Table II). 
The eccentric stations were within 2 m of the main station 
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and marked by masonry nails driven into the concrete. The 
geodetic position of each eccentric station was computed on 
a HP 9815 calculator using National Geodetic Survey (NGS) 
geodetic programs. Initial pointings for the T-2 
theodolites were to MUSSEL from stations SQUARE, CONK, and 
USE MON; and to MONTEREY RADIO STATION KMBY MAST from 
station GEOCEIVER . 







TABLE 


II 






GEODETIC CONTROL POSITIONS 




Station 


X (m) 


Y (m) 


Latitude 


Longitude 


SQUARE 
SQUARE ecc 


7974.86 

7976.10 


3909.43 

3907.89 


36°37 ' 07 . 175” 
36°37’07. 140" 


121°51 ' 00 .276" 
121°51 ' 00 . 230" 


CONK 

CONK ecc 


6978.19 

6976.70 


2828.77 

2827.53 


36°36' 32. 130" 
36°36 ' 32 . 100" 


121°51 ' 40 . 397" 
121°51 ' 40 . 460" 


USE MON 


5598.98 


1-982.76 


36° 36 ' 04 . 685" 


121°52' 35.900" 


GEOCEIVER 
GEOCVR ecc 


4371.50 

4372.99 


2840.28 

2839.97 


36 ° 36 ' 32.512" 
36°36 ' 32 . 510" 


121°53 ' 25.286" 
121°53 ' 25 . 230" 


MUSSEL 


3220.17 


4247.23 


36°37 ' 18 . 151" 


121°54' 11.628" 


KMBY MAST 


3641.23 


3588.23 


36°36 '56.789" 


121°53 ' 54 . 678" 



The data acquisition procedures were similar to normal 
survey operations for a large-scale hydrographic survey 
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using visual control with the addition that on each fix four 
ranges were recorded automatically by the AUTOCARTA II 
system. The four theodolite observers trained their 
instruments on the vessel as it moved along the track line. 
Position fix marks were communicated over voice radio and 
the T-2 directions to the vessel were manually recorded by 
the observers. 

To eliminate systematic errors associated with the Del 
Norte system, the Trisponder units were calibrated over a 
measured geodetic base line 2565.897 m long (CONK 1984 to 
REY 1984) on 24 and 25 November 1984. A base-line 
calibration prior to survey operations is a standard field 
procedure for this type of equipment. Daily dynamic system 
checks were made similar to NOS procedures [Holder, 1983]. 
The observed "drifts” in the Trisponder units were recorded 
from an ending base-line calibration on 7 December 1984 and 
were applied as range correctors during the data processing 
phase of the experiment. 

Weather conditions during the 2 days in which horizontal 
directions to the vessel were recorded was fair. The 
theodolite observers had difficulty in maintaining an 
adequate initial pointing due to winds gusting to 30 knots 
which caused vibration of the tripods. Initial checks were 
made at the end of each track line and recorded. Observed 
changes in the initial directions were recorded and applied 
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during the data processing phase to eliminate any systematic 
errors in the angles. On the first 2 days of field 
operations, a total of 176 positions, each with four ranges 
and four azimuths, were recorded. An additional 175 
positions with four ranges were acquired using the AUTOCARTA 
II system during the last 3 days of the field work. 

B. DATA ACQUISITION SYSTEM 

The AUTOCARTA II (Figure 3.2) is a portable automated 
hydrographic data acquisition system. Desired vessel tracks 
can be preprogrammed into the computer to alleviate the 
necessity of steering the vessel along range arcs for track 
control. A real-time graphic copy of the vessel’s actual 
course was plotted from Autocarta II output on a Houston 
Instruments Model DP- 3 Plotter. A Model 540 DDMU recorded 
the four ranges at 1-minute intervals along the track line. 
Input and output • communication to the AUTOCARTA II were via 
a modified TI-743KSR keyboard/printer. 

C. PERFORMANCE OF DATA ACQUISITION SYSTEM 

The operational difficulties experienced during the data 
acquisition phase of the field work could be described as 
"standard” for a hydrographic survey of this type. To 
ensure adequate backup was available for Trisponder units, a 
total of seven remote and two master units were on hand. 

Two of the remote units did not responded when interrogated 
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Figure 3.2 Data Acquisition System. 







during the base-line calibration of the system and had to be 
removed for the entire study. New deep- cycle marine 
batteries were used as power to the 24-volt Trisponder 
system. These batteries were charged each night to ensure 
sufficient power during the survey. 

A cable linking the AUTOCARTA II to the keyboard/printer 
was not shipped with the equipment. Therefore, the data 
were acquired in the normal non- automated fashion by hand 
recording the range rates as they were displayed on the DDMU 
on the first day of the survey. The needed cable was 
available by the next day and the complete automated system 
was operational for the rest of the survey. 

The AUTOCARTA II was set up to interrogate the remote 
units on shore from left to right when facing shoreward from 
the vessel which is standard NOS convention. On the first 
day of data acquisition, day 332, very erratic range rates 
were being observed on the DDMU from station USE MON 1978 
(remote code #76). The update rate was set at 1-second 
intervals and over half of the time either a totally 
ambiguous rate or no rate at all was recorded on the DDMU. 

In normal field operations , a remedy to the situation is to 
raise or lower the remote unit to eliminate the possibility 
of "null zones" and "skip zones.” The remote unit on 
station USE MON 1978 was raised and lowered with no change 
in the performance of the Trisponder observed. 
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IV. MATHEMATICAL MODEL 



A. LEAST SQUARES MODEL 

To determine the ’’best" estimate of a position, 
redundant observations are necessary so that the method of 
least squares can be used. The least squares adjustment by 
observation equations (variation of coordinates) was 
selected for this thesis (rather than least squares 
adjustment by condition equations) for reasons of simplicity 
and clarity [Bomford, 1980, pp . 127-128]. 

A basic discussion of the method of least squares by 
observation equations is given below [Bomford, 1980; Mikhail 
and Gracie, 1981; Cross, 1981; Heinzen, 1977; Kaplan, 1980; 
Saxena, 1972; Silva, 1982]. 

1 . Unweighted Equations 

The equation derivations will be initially 
accomplished for measurements of equal accuracy. 

Let L^ represent independent observed quantities; v^ 
be residuals of the observed quantities; and x, y, and z be 
the unknown coordinates of a point in a rectangular 
coordinate system. Then for each observation, i, an 
observation equation is formed such that 

L i + v i = f i (x,y,z) (4.1) 
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where is a linear or non-linear function. When is a 
non-linear function, a Taylor series is formed about the 
approximate values of x Q , y Q and z Q such that the second 
and higher order terms are small enough to be neglected. 
The observation equations then become 



If - f ^ (x Q ,y , z Q ) , which is the difference between 
the observed quantities and the computed approximate values 
of the unknowns , then 




f i ( x 0 +dx » y 0 +d y> z o +dz > • • •) 



(4.2) 



or 



L i + Vf = f i (x Q +dx, y Q + dy, z Q +dz) 

+ a^dx + b^dy + c^dz + . . 



(4.3) 



v^ = a^dx + b^dy + c^dz + . . • ~ P± 



(4.4) 



where 



x = x Q + dx 



(4.5) 



y = y Q + d y 



(4.6) 



z = z Q + dz 



(4.7) 



and 



= 3f/3x, b i = 3f/3y> = 3f/3z 



(4.8) 
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so that in matrix form the observation equation 4.4 becomes 



V = AX - L (4.9) 

which ultimately yields the normal equations 

NX - U = 0 (4.10) 

where 

N = A T A (4.11) 

and 

U = A T L (4.12) 



Then substituting into Equation 4.10 and rearranging terms 
gives 

X = N _1 U = (A T A) -1 A T L (4.13) 

which is valid for observations of equal weight. 

In actual surveying, computed observations differ 
from actual observations by a value known as a residual, v. 
The residual is the difference between an observed value, 
x..- and the estimate of the true value, x, so that 

x = Xj_ + Vj_ (4.14) 

The basic condition of the least squares method is 
that the sum of the squares of the residuals equal a 
minimum. Therefore, 
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minimum 



(4.15) 



£ (vi ) 2 



where 




. + v 



n 



2 



for values of equal weight. 



2 . Weighted Equations 



When observations are made utilizing various methods 



and systems, weights must be applied to account for the 



more precise observations. The relationship between weights 



where k is a constant of proportionality having the same 
value for all observations in one system when the different 
standard errors of the observations are uncorrelated and of 
the same precision, as is the case for range-range, 
azimuth-azimuth, and range-azimuth positioning systems. 

Thus 



p 

lower standard error, <j 5 and lower variance, <j , of the 




w = k / ct^ 



(4.16) 



k = 1 



(4.17) 



Then 




(4.18) 
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which in matrix form is the diagonal matrix 

2 



W = 



l/Gl 



l/(7 2 ‘ 



1/a 



n 



For weighted observations, Equation 4.15 becomes 



(4.19) 



£ w^(v^) = minimum 



(4.20) 



where 

£ W^(v^^) = WjV^ + W£V£ 



2 

Vo + . . . + W^V 



1 1 



or in matrix form 



T 

V WV .= minimum 



(4.21) 



which after substitution and differentiation with respect to 
X becomes 



X = (a t wa)' 1 a t wl 



(4.22) 



This system of equations is then solved for the 
matrix X elements which are then added to a provisional set 
of coordinates at a point, yielding a new pair of 
coordinates. This process is iterated until the values of 
the X elements become small enough to meet a specific 
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tolerance. At this point the method is said to have 
converged. What remains is the "best” estimate of x and v, 
designated as x and v. 



B. PRECISION OF ADJUSTED VALUES 

The first right-hand term of Equation 4.22 is the 
variance- covariance matrix, 

(A T WA) ~ 1 = N" 1 = Q (4.23) 

and it can be shown that , 



Q = (A T WA) * 1 = 

The variance- covariance matrix of the least squares estimate 
of x, cr(&) , is 



^11 ^12 
^21 c *22 



• • 



e © t • 



« o • • 



l nn 



(4.24) 



a(x) = a 0 2 (A t WA) 1 = a Q 2 N' 1 = a 0 2 Q 



a(x) = 



°x CT xy 



CTyx 



'n 



(4.25) 
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is known as the variance of an observation of unit 



p 

Here a 0 



weight represented as 



CT 0 2 = V T WV / (n - m) 



( 4 . 26 ) 



The standard error of an observation of unit weight, a-, is 

U > 

therefore 



where n is the number - of observations, and m is the number 
of unknowns. The term (n - m) is known as the degree of 
freedom. 

Consequently, the standard error of the adjusted 

r 

coordinates, <j x and <jy, are given by; 




( 4 . 27 ) 



CT o V ^xx 



( 4 . 28 ) 




( 4 . 29 ) 



and the standard position error is 




( 4 . 30 ) 
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Finally, the semi-major (a a ) , semi-minor (at>)> and 
orientation (0) axes of the error ellipse are given as 



°a a o \J Qmax (4.31) 

<*b = a o 4 Qmin • (4.32) 

where 

^raax.min “ (-*-/^) (^xx + Qyy) (4.33) 

± (1/2) y(Qxx-Qyy) 2 + 4 Qxy 2 

and 

6 = (1/2) arc tan [ 2 Q xy / (Q xx - Q yy ) ] (4.34) 
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V. DATA ANALYSIS AND RESULTS 



A. DATA PROCESSING PROCEDURES 

Data acquired in the field were manually logged into 
data files on the IBM 370/3033AP mainframe computer at NPS 
via .the VM terminals. These data were processed using 
computer programs (Appendix A) written to be compatible with 
the VS FORTRAN compiler and used many of the existing 
library subroutines. The main program which computes the 
least squares position from three- or four-range LOPs , 

SLEAST (Appendix B), was composed by Lt . Nick Perugini, 

NOAA, and the author. The article by Cross [1981] covers 
the subject in detail. The least squares programs for the 
azimuth and hybrid MLOP methods are FORTRAN programs LSQAZ4 , 
R3+AZ, R2+AZ2, and RR+AZ (Appendix A). These programs are 
modifications of programs SILVA1 and SILVA3 , written by 
Silva [1982] . Computation of coordinates and accuracies 
from two LOPs were made by modifying the program UCOMPS 
(Appendix A) into RR2L0P , AZ2L0P , RAZ2L0P . 

Processing procedures commenced by correcting the 
observed ranges and azimuths for previously determined 
systematic errors prior to logging the data into manageable 
files. The least squares computer programs output the 
position of the vessel (in X and Y coordinates based on a 
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Modified Transverse Mercator (MTM) projection^); the final 
residuals; and the computed accuracy, Cp> of each position 
fix. The output from the processing programs for two LOPs 
were the position coordinates of the vessel, the angle of 
intersection of the LOPs, p, and the standard error of 
position, CTp , as determined in Chapter II for each position 
method . 

The geodetic positions of the shore stations used for 
hydrographic control were converted to MTM coordinates via 
the utility program UCOMPS (Appendix A) . 

B. DATA ANALYSIS 

The main objective of the thesis was to investigate 
whether positions computed from MLOP were more accurate than 
positions determined by only two LOPs. Models and 
algorithms were developed so field data could be 
appropriately analyzed. For each positional fix, eight LOPs 
(four ranges and four azimuths) were measured 
simultaneously. The eight LOPs were processed separately 
into ten different sets of positioning methods. Depending 



^The MTM projection is a Universal Transverse Mercator 
projection except the Central Meridian is near the survey 
area instead of being at a preselected meridian [Wallace, 
1971, pp . 3-4]. A Central Meridian of longitude 121 o 53'00 M W 
was chosen for this thesis. The X coordinate assigned to 
the Central Meridian, the FEST, was chosen as 5,000 m. The 
Controlling Latitude, CLAT , or the distance in meters from 
the equator to some reference latitude, was computed to be 
4,050,000 m. 
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on the method, some of the ten sets were separated into 
various subsets by using different combinations of 
observations from the four shore stations. 

Let the four shore stations from which the LOPs 
originated be designated A, B, C, and D. If four-range (or 
four-azimuth) LOPs are to be used for positioning, there is 
only one unique combination of four-range (or four-azimuth) 
LOPs which can solve the problem. Using only three LOPs, 
there exists four possible combinations of range (or 
azimuth) LOPS--ABC, ABD , ACD, and BCD. 

The problem becomes more confusing with hybrid 
combinations of range LOPs and azimuth LOPs. A practical 
limitation was imposed by deciding that each azimuth LOP 
must originate from the same station (or eccentric station) 
as the range LOP. Thus, if three LOPs were to be determined 
from two ranges and one azimuth, the azimuth observation 
must be made at one of the stations (or its eccentric 
station) where the ranging unit is located. For example, if 
Trisponders are set on stations A and B (or the eccentric 
stations AE and BE). The azimuth LOP must originate from 
either station A or station B. There are a total of 12 LOP 
combinations which must be considered when using three 
ranges and one azimuth, or two ranges and one azimuth 
(Table III). 
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LOP 


TABLE 

COMBINATIONS FOR 


III 

TWO HYBRID METHODS 


3 Ranges 


& 1 Azimuth 


2 


Ranges & 


1 Azimuth 


3 Ranges 


Azimuth 


2 


Ranges 


Azimuth 


ABC 


A 




AB 


A 


ABC 


B 




AB 


B 


ABC 


C 




AC 


A 


ABD 


A 




AC 


C 


ABD 


B 




AD 


A 


ABD 


D 




AD 


D 


ACD 


A 




BC 


B 


ACD 


C 




BC 


C 


ACD 


D 




BD 


B 


BCD 


B 




BD 


D 


BCD 


C 




CD 


C 


BCD 


D 




CD 


D 



All possible combinations of the 10 data sets were 
processed resulting in a total of 56 positions being 
computed for most of the fixes (Table IV). 

Gaps exist between consecutive fix numbers in the output 
files (Appendices B, C, and F) due to spurious ranges or 
unresolved blunders which were rejected. Data were rejected 
if the standard error of position exceeded 3<j for that LOP. 
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TABLE IV 

DATA FILE COMBINATIONS 



Method Possible Combinations 



4 LOPs 


4 R 




i 


4 Az 




i 


3 R & 1 


Az 


12 


2 R & 2 


Az 


6 


3 LOPs 


3 R 


- 


4 


3 Az 




4 


2 R & 1 


Az 


12 


2 LOPs 


2 R 




6 


2 Az 




6 


1 R & 1 


Az 


.4 


TOTAL 
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The following assumptions were made concerning the 
overall accuracy determination of each position: 

(a) All LOPs were independant and uncorrelated. 

(b) All data were free of blunders and systematic errors. 

(c) The standard error of range LOPs was 3 m at all 
ranges . 

(d) The standard error of azimuth LOPs was 4 m. 



61 



(e) The measured ranges were referred to the same point 
of time during interrogation by the DDMU time-deskew 
function [Anonymous, 1982]. 

(f) Measurements by the azimuth observers were all of 
equal weight . 

1 . Range Measurement Methods 

An analysis was made of the acquired range data by 
the number of LOPs . The mean value of <7 from each 

F 

positioning method is illustrated in Table V. 

a. MLOP 

FORTRAN program SLEAST (Appendix B) processed 
the three- and four-range LOPs (Appendix C) using the least 
squares method. The mean <7 value for each day of data and 

F 

the minimum and maximum <7 were computed (Table V). For 

F 

the 3 days in which range measurements were acquired, an 
average <7 of 1.3 m was attained with the four LOPs and 1.5 

F 

m with the three LOPs. The geometric accuracy of the MLOP 
fixes exceeded the standard error of the Del Norte system by 
at least 1.5 m. The slight variation in the mean <7 ’ s for 

F 

the three LOPs was due to the geometrical variation of the 
ranging net from each shore station. 

b. Two LOPs 

The two-range LOPs were processed by program 
RR2L0P (Appendix A) to compute the coordinates of the vessel 
plus the standard error of position, <7 , for each fix. 

F 

Analysis of the data shows a direct relationship between the 
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TABLE V 

POSITIONAL ACCURACIES FOR RANGE LOPS 






dp 


CTp 


CTp 


ANGLE OF 


LOPs 


DAY 


MIN 


MAX 


MEAN 


INTERSECTION 






(m) 


(m) 


(m) 


(2 LOPs) 


4 R 


332 


0.3 


00 

C\1 


1.2 




4 R 


333 


0.4 


2.9 


1.4 




4 R 


334 


0.5 


2.4 


1.3 




3 R 


332 


0.1 


3.6 


1.1 


"V 






0.04 


4.5 


1.7 




3 R 


333 


0.04 


5.7 


1.3 








0.01 


3.7 


2.0 




3 R 


334 


0.1 


2.4 


0.8 








0.4 


4.5 


2.2 


*\ 


2 R 


332 


4.8 


7.9 


6.0 


147° 






7.7 


236.7 


29.2 


2° 


2 R 


333 


4.4 


8 . 1 


5.7 


28° 






4.4 


578.2 


12.3 


1° 


* NOTE: 


For positions 


which 


could be 


computed by more 


than one 


combination of 


LOPs 


( ® • S • » 4 


combinations of 


3 -range 


LOP fixes exist 


when 


4 ranges 


have been 


measured) , the 


upper line for 


each day 


shows the 


combination of 


LOPs giving the minimum 


mean <jp ' s 


and the 


lower line for 


each day shows 


the combination 


of LOPs 


giving 


the maximum mean cTp ' s . 
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angle of intersection of the LOPs and cjp , as would be 
expected. A wide variation in the computed <Jp ’ s exists due 
to the varying angle of intersections. The pre-selected 
sounding line orientation was programmed into the AUTOCARTA 
II prior to the acquisition of data. Since the purpose of 
the field work was to determine the position of the vessel 
using MLOP , the angle of intersection of each individual 
pair of LOPs was ignored. Thus at fix number 19 the 
computed <jp was 237 m due to an angle of intersection of 
1.6°. During normal survey operations this situation would 
not be accepted and alternate plans would be made to obtain 
coverage in that area. Herein lies a benefit of using MLOP 
since the position of the vessel can be determined from 
selected shore stations which yield' the best geometric 
configuration. As could be expected, the required accuracy 
of at least 5 m was not met when using two range LOPs with a 
standard error of 3 m (Table V). 

Overall the capability of acquiring MLOP 
automatically from measured ranges has proven to be a very 
viable system. Mean accuracies from 1.2 to 1.6 m were 
observed depending on the number and geometric combinations 
utilized. With these accuracies all nearshore positioning 
requirements can be met. 

Using range MLOP for hydrographic positioning 
has many advantages. The capability of running a fully 
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automated survey with systems such as the AUTOCARTA II not 
only increases the productivity but also reduces blunders , 
thereby increasing the quality of the survey. Once remote 
transponders are erected over the control stations the need 
for attendant operators is negated thus reducing the number 
of personnel required to acquire data. The portability of 
the system (except possibly for the batteries) is an 
additional advantage in rugged terrain. 

2 . Azimuth Measurement Methods 

The use of an angle measuring instrument , such as a 
T-2 theodolite, for positional control of sounding lines is 
a laborious task. Continuous tracking of the vessel is 
difficult because the horizontal tangent screw has a finite 
range. Observers can be tracking the vessel with the 
tangent screw only to find the end of the drive mechanism 
seconds before the position fix. Thus the problem of 
tracking a vessel who’s aspect is changing rapidly with 
respect to a nearby observer is magnified. Older T-2 
theodolites have an inverted image of the target which may 
confuse an observer since the vessel seems to be moving in 
the opposite direction. Also the reading of the observed 
direction is through a secondary scope. Therefore, an 
observer must stop tracking the vessel, align the horizontal 
plates, read the direction and record the value. The lack 
of digital telemetering capabilities requires the observed 
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value to be relayed to the survey vessel by voice radio. 
Obviously, there exists a potential for many sources of 
errors if position fixes are required at frequent 
(30-second) intervals. 

Using nonautomated systems (i.e., visual azimuths) 
restricts position fixing of the vessel to well spaced 
intervals (minimum of 30 seconds). Intermediate positions 
along the sounding line must be interpolated via 
dead- reckoning . In most cases, the actual track of the 
survey vessel does not exactly correspond to the plotted 
survey line. On a 1 : 5 , 000- scale survey, this additional 
error could exceed the specified accuracy requirements. 

Although the above limitations seem to make the 
method too complicated and error prone for hydrographic 
control, it is still one of the most exact positioning 
methods available for large-scale surveys. Experienced 
observers have little difficulty in acquiring accurate data. 

a. MLOP 

Even though a 4-m standard error was applied to 
all azimuth LOPs , the least squares position from three or 
four azimuths proved to be the most accurate method 
investigated. Examination of the tabulated <jp ' s (Table 
VI) shows a mean value of 0.7 m on day 332 and 1.0 m on day 
333 when the four LOPs were used. The mean <jp ' s from three 
LOPs ranged from 0.7 m to 1.3 m over all combinations from 
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TABLE VI 

POSITIONAL ACCURACIES FOR AZIMUTH LOPS 







a P 


CTp 


a P 


ANGLE OF 


LOPs 


DAY. 


MIN 


MAX 


MEAN 


INTERSECTION 






(m) 


(m) 


(m) 


(2 LOPs) 


4 Az 


332 


0.3 


1.7 


0.7 




4 Az 


333 


0.1 


2.9 


1.0 




3 Az 


332 


0.00 


2.5 


0.7 








0.01 


9.5 


1.2 




3 Az 


333 


0.00 


4.9 


0.8 


«IU 






0.05 


4.3 


1.3 




2 Az 


332 


5.7 


8.9 


6.3 


41° 






6.3 


423.6 


26.7 


1° 


2 Az 


333 


5.7 


12.9 


6.4 


26° 






5.7 


84.2 


8 . 6 


4° 


* NOTE: 


For positions 


which 


could be 


computed by more 


than one 


combination o 


f LOPs 


(e.g. , 4 


combinations of 



3-azimuth LOP fixes exist when 4 azimuths have been 
observed), the upper line for each day shows the 
combination of LOPs giving the minimum mean <jp ' s 
and the lower line for each day shows the combination 
of LOPs giving the maximum mean <jp ' s . 
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the 2 days of data. The accuracy capabilities of the method 
are evident . 

b. Two LOPs 

Although the least squares adjusted positions of 
the azimuth MLOP proved to be extremely accurate, the 
two-azimuth-LOPs condition suffers from the same degradation 
in accuracy as the two - range-LOPs case. The error in 
position, ct s in cm, was determined from Equation 2.46. 

Since the standard error of an azimuth LOP has been defined 
as 4 m throughout this thesis, must be evaluated. From 
the standard conversion formula (Equation 2.18) then 

<j a = ( b / R ) 206265 s/rad (5.1) 

where R is the distance in m, and b is the standard error of 
4 m. Since the azimuth LOPs result in Xp and Yp , the 
coordinates of the vessel, then the distance from the 
computed coordinate to each shore control station is 

dl = J ( X L - Xp ) 2 + (Y r Yp) 2 (5.2) 

and 

d2 = J (X 2 -Xp) 2 ♦ (Y 2 -Yp) 2 (5.3) 

Rather than compute separate error values for each distance 
an average of the distance from both stations, d a , was used 
where 

d a = (dl + d2) / 2 (5.4) 
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and Equation 5.1 becomes 



<J a = ( 4 / d a ) 206265 (5.5) 

which can be substituted into Equation 2.46 to compute a s 
(in cm). Consequently, the computed <jp (<j s in cm times 
100) is a function of both the range to the vessel and the 
angle of intersection of the two LOPs. 

Analysis of all of the (Tp ' s attained in this 
fashion shows a mean value of about 6.3 m for all data 
acquired (Table VI). Although this is rather high for an 
error in position from two accurate azimuths, the cause can 
be easily explained. The initial assumption that each 
azimuth LOP was in error by 4 m at all ranges is rather 
excessive for theodolite measurements. A more realistic 
investigation of the angular error of dynamic azimuth 
measurements should be considered in the future. Standard 
field procedure for azimuth hydrographic control is to 
ignore the stated random error of the theodolite because the 
measured directions are only read to the nearest minute of 
arc. The only errors associated with a position determined 
by azimuths are the initial correctors and any blunders 
discovered during data acquisition and post processing. 

3 . Hybrid Methods 

Numerous combinations of control configurations 
exist when ranges and azimuths are mixed together. The data 
acquired were processed by FORTRAN programs R3+AZ, R2+AZ2, 
and RR+AZ (Appendix A) . 
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a. Three-Range and One- Azimuth LOPs 
Three-range and one- azimuth LOPs exhibit larger 

mean values of (jp ' s than in the prior MLOP methods (Table 
VII). One explanation for this difference is the geometry 
of the fix has a greater effect on the hybrid methods than 
on the other methods of position control. In this 
situation, one intersection is at 90° and up to four 
intersections are at arbitrary angles depending on the 
position of the vessel relative to the shore stations. If 
the angle of intersection of the range LOPs is too large 
(greater than 150°) or too small (less than 30°) then the 
quality of- the fix may be degraded. The 11.8-m and 12.5-m 
maximum values for the C 7 p ’ s observed on both days (days 332 
and 333, Table VII) were found using configuration ABDD 

r 

(three ranges from stations A, B, and D and the azimuth from 
station D) . In both instances, the angle of intersection 
between the LOPs from stations A and B was greater than 
150°. Therefore, the hybrid method is very sensitive to the 
angle of intersection of each LOP pair where a range and an 
azimuth are not taken from the same station. 

b. Two-Range and Two-Azimuth LOPs 

Analysis of the results obtained from processing 

the two -range and two -azimuth LOPs found not as much 

variability in the range of the CTp ' s as in the (jp ’ s of 

three-range and one-azimuth LOPs. An overall mean ov, of 

F 
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TABLE VII 

POSITIONAL ACCURACIES FOR HYBRID LOPS 






a P 


a P 


<7p 


ANGLE OF 


LOPs 




DAY MIN 


MAX 


MEAN 


INTERSECTION 






(m) 


(m) 


(m) 


(2 LOPs) 


3 R & 1 


Az 


332 0.1 


11.8 


1.1 


-A. 






1.2 


6.3 


3.0 


/> 


3 R & 1 


Az 


333 0.1 


6.7 


1.7 








1.2 


12.5 


5.5 




2 R & 2 


Az 


332 0.3 


2.2 


0.8 








1.2 


5.3 


2.1 




2 R & 2 


Az 


333 0.2 


3.5 


1.4 








0.2 


5.0 


2.4 




2 R & 1 


Az 


332 0.01 


5.1 


0.8 








0.3 


4.9 


1.9 




2 R & 1 


Az 


333 0.02 


4.2 


1.1 








0.1 


7.7 


2.6 


** 


1 R & 1 


Az 


332 3.0 


5.0 


5.0 


90° 


1 R & 1 


Az 


333 5.0 


5.0 


5.0 


90° 


* NOTE: 


For positions 


which could be 


computed by more 


than one combination of 


LOPs (e 


•g* > 4 


combinations of 


3 - range 


LOP 


fixes exist 


when 4 


ranges 


have been 


measured ) , 


the upper line for each day 


shows the 


combination 


of LOPs giving the i 


minimum 


mean <jp ' s 


and the 


lower line for 


each day 


shows 


the combination 


of LOPs 


giving the maximum mean 


CTp s . 
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1.5 m on day 332 and 1.9 m on day 333 was observed. 

Observing both a range and an azimuth, which intersect at 

90°, at each of the two stations causes the cr values to be 

F 

more consistent. When the two normal LOPs from one station 
are combined with two normal LOPs from another station a 
strong position is determined. Since an angle of 
intersection of 90° is the best possible geometry for a 
two-LOP fix, then this position is made from two accurate 
sets of LOPs. 

c. Two-Range and One- Azimuth LOPs 

Overall mean values of <j of 1.4 to 2.0 m were 
determined when two-range and one-azimuth LOPs were used for 
positioning. Variability exists in the range of Qp values 
observed due to the geometry problems stated previously. 

Two of the LOPs intersect at 90° and then a third LOP is 
added to this accurate configuration to determine the 
three-LOPs fix. 

d. One-Range and One- Azimuth LOPs 

The observed <j for the two-LOP Range and 
Azimuth method varied very little due to the minute change 
in the intersection angle of the LOPs. Since the ranges 
were measured from eccentric stations for three of the four 
control sites, the angle of intersection was slightly 
different than 90°. A value of 5 m for the standard error 
of position was determined for most of the position fixes. 
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Due to the consistency of the intersection angle 
of the LOPs, the range-azimuth LOPs method is prefered for 
large-scale hydrography. With the present interpretation of 
the IHO specifications, the 5-m cjp would be acceptable for 
positional accuracy. 

C. COST COMPARISONS OF POSITIONING METHODS 

The increased cost of using MLOP may outweigh the 
benefits of an increase in accuracy. Conventional two-LOP 
methods require the minimum equipment and personnel to 
produce a hydrographic survey. The continuing increase in 
marine positioning requirements may force the hydrographer 
of the future to decide what price should be paid to produce 
a more accurate survey. 

Component costs of the equipment used in this thesis 
were obtained in June 1985 from manufacturer representatives 
(Table VIII). Prices quoted are the government contract 
rate in U.S. dollars. Table IX combines the costs of all 
the components depending on the method used to position the 
vessel. Only the four-LOP and two-LOP methods -were 
addressed as the costs of the other methods can be inferred. 
The costs for the two-LOP method assumes the data will be 
acquired by the conventional nonautomated procedure commonly 
used on large-scale surveys. Equipment which is common to 
all methods (i.e., vessels, tripods, signals, etc.) have 
been ignored in the tabulation. 
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TABLE VIII 



COMPONENT EQUIPMENT COSTS 



COMPONENT COST ( $ ) 

Del Norte Trisponder: 

Model 542 DDMU 14 , 500 

Master 9 , 500 

Remote 9,600 

RACAL-DECCA AUTOCARTA II 24,750 

Houston Instrument Plotter 6,500 

Wild T-2 Theodolite 6,500 

Portable VHF Radio 800 



The price of equipment for the four-range ML0P automated 
method approaches $100,000. In comparison, equipment for 
the four- azimuth ML0P method is about one- third the price of 
the four- range ML0P method. The difference in equipment 
costs in the two methods is very deceptive. To acquire four 
azimuths, four observers are needed to track the vessel. 

The cost for this additional personnel over the automated 
ranging system could approach the difference in equipment 
costs in just over one year. The labor costs involved would 
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TABLE IX 

COST COMPARISON OF EQUIPMENT 



Four-Range 


MLOP 




540 DDMU 
Mci s ter 

Remotes (4 (3 $9,600) 
AUTOCARTA II 
Plotter 




$14,500 
9 500 
38,400 
24,750 
6,500 




TOTAL 


$93,650 


Four-Azimuth MLOP 




T-2 Theodolites (40$ 
Portable Radios (4 (a $ 


6.500) 

860 ) 


26,000 

3,200 




TOTAL 


$29,200 



Two-Range & Two-Azimuth MLOP 



540 DDMU 
Master 

Remotes (2 0 $9,600) 
AUTOCARTA II 
T-2 Theodolites (20$ 
Portable Radios (20$ 


6,500) 

860 ) 


14,500 
9,500 
19 200 
24 7 50 
13,000 
1*600 




TOTAL 


$82,550 


Two -Range 


LOPs 




540 DDMU 
Master 

Remotes (2 0 $9,600) 
AUTOCARTA II 




14,500 

9,500 

19,200 

24^750 




TOTAL 


$67,950 


Two-Azimuth LOPs 




T-2 Theodolites (2 0$ 
Portable Radios (20$ 


!6,500) 

;86o) 


13,000 

11600 




TOTAL 


$14,600 


Range and Azimuth LOPs 




540 DDMU 

Master 

Remote 

T-2 Theodolite 
Portable Radio 




14,500 
9,500 
9 600 
6 500 
800 




TOTAL 


$40,900 
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depend upon the skill and experience of the personnel used 
to acquire the data. For four observers paid $8.00 per 
hour, for an 8-hour day, for 260 days per year, the increase 
due to labor would be $66,560 (4 x 8 x $8 x 260) per year. 
Added to the equipment costs of $29,200, the total price 
exceeds $95,000. 

Another expense associated with the four-azimuth MLOP 
method would be the cost of the hardware needed to process 
the data using a least squares algorithm. A data processing 
computer would be necessary to perform the repeated 
iterations required by the algorithm. Additionally, a major 
disadvantage of the nonautomated data acquisition method is 
that observed azimuths must be hand logged into the data 
processing system. This manual intervention increases the 
existance of numerous blunders and involves additional labor 
costs which could approach the cost of acquiring the 
four- azimuth LOPs. 

With these facts in mind, the overall costs of the 
four-azimuth method could approach the cost of the automated 
four-range method. Given this circumstance, the optimal 
system would be to use an automated ranging method to 
increase production, reduce blunders in recording the 
observations , and have better control of the survey 
operations . 
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D. OPTIMIZED SYSTEM 



One of the objectives of this thesis was to determine an 
optimal system for the acquisition of positional data in 
hydrography. Results from the data analysis proved 
conclusively that MLOP do improved fix accuracy over two 
LOPs . Every effort was made to duplicate the actual field 
work currently done by NOS so that a valid comparison of 
methods could be made. 

To adequately discuss an optimal system a few 
considerations must be made. First, a consideration of the 
accuracy requirements must be made. This accuracy 
requirement will be dependent upon the purpose of the 
survey. Even with the initial accuracy estimates used for 
this thesis, all of the MLOP and the range-azimuth two-LOP 
methods will meet the 5-m positional accuracy. Therefore, 
the optimal system would be dependent upon the final product 
the hydrographer wishes to produce. If automation is the 
overriding criterion, then a ranging MLOP method would 
adequately fulfill the needs of the survey. An automated 
system, such as AUTOCARTA II, can be bought or leased 
off-the-shelf, to provide the user with an automated MLOP 
system. Automation frees personnel to do other tasks since 
the need for shore observers is eliminated (except for 
station maintenance). By acquiring the data automatically 
on board the survey vessel, the existance of blunders are 
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minimized, thereby decreasing the processing time. Since 
the capability exists for pre-selecting track' lines and 
controlling the movements of the vessel, more production is 
possible with an automated system. Overall the use of range 
MLOP and an automated DAS is the most feasible method to 
conduct accurate near-shore hydrography. 

Second, if the purpose of the survey dictated that the 
accuracy of position must be to within 1 m, or less, then 
more effort must be made to acquire this accuracy. Even 
when a 4-m standard error was used for each LOP, this thesis 
found that the <j for 63% of the positions determined by 

r 

four- and three- azimuth LOPs were less than or equal to 1.0 
m. Although the method is susceptible to numerous errors, 
is labor intensive, and logistically complex, one could 
produce a survey of a nearshore area to within 1.0-m 
accuracy with good field procedures. 

Third, cost must be considered. The benefit of using 
just two LOPs is that minimal equipment and labor costs are 
involved to acquire adequate data. The use of MLOP will 
increase the initial cost of equipment at least twofold. 

The cost may be prorated over time so in 4 to 5 years the 
benefit of producing a higher quality product would 
compensate the hardware costs. Additionally, the cost to 
establish more shore control must be considered. More LOPs 
would mean more control stations are required than presently 
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needed for conventional hydrographic surveys. Detailed 
survey planning could minimize the cost by choosing the 
optimum number of stations needed for the method of data 
acquisition . 
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VI. CONCLUSIONS 



A. MLOP FOR ON-LINE CALIBRATION OF RANGING SYSTEMS 

Short-range microwave survey systems, such as the Del 
Norte Trisponder, require a static base-line calibration 
over a known geodetic distance to eliminate systematic 
errors. Daily system checks must also be made during the 
course of the survey to detect any calibration drifts over 
time. Daily system checks can be made in either a static or 
dynamic mode. Static checks can be made by taking fixed 
point observations alongside a known point, three-point 
sextant fixes with a check angle, theodolite intersections 
from shore control, or range-azimuth positions using a total 
station EDM unit. Dynamic checks can include steering a 
range with predetermined crossing angles or underway 
simultaneous theodolite intersections [Holder, 1983]. 

Experience has shown that whatever system is used to 
perform daily calibration checks of the equipment, a 
substantial amount of time is consumed on this exercise in 
quality control. Weather permitting, daily system checks 
are made prior to and after each day’s hydrography. Under 
ideal conditions, the amount of time consumed to obtain the 
required calibrations usually is 1 hour of the survey day. 
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A distinct advantage of using MLOP for position control 
in hydrography is that redundant data are acquired 
simultaneously on each fix. By having a third or fourth 
LOP, errors in position could be detected more readily than 
in the conventional survey mode where everything was assumed 
to be correct until it was too late. Thus the need to halt 
survey operations to make a daily system check could be 
eliminated. Suppose four-range MLOP were to be used to 
position the vessel during a survey. Once the system was 
on-line and warmed-up, a number of positions could be 
obtained by the automatic system while underway to the 
working grounds. If a large standard error of position, 

(jp , was determined by the software within the computer, 
then an examination of the residuals for each LOP could be 
examined to determine which range or ranges were in error. 
This error could be from equipment malfunctions, errors in 
the locating of shore stations , or erroneous ranges due to 
low signal strength, reflections, etc. Additional fixes 
could confirm, or reject the existance of a substantial 
error. Once isolated, the unit (or units) causing the error 
can be removed or replaced. Here the advantages of having 
four-range MLOP as control can be realized. Should one of 
the units fail, there would still be a redundant range 
acquired allowing for a least squares computation. If for 
some reason two units failed, operations could continue 
using conventional two-range LOPs . 
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B. FUTURE ASPECTS AND CONSIDERATIONS 



The use’ of MLOP for hydrographic sounding line control 
should be the primary method for future nearshore surveys . 

As advances in electronic technology continue, the economy 
of miniature computer installations will make automated 
hydrographic data acquisition systems more feasible. NOS is 
presently in the process of revitalizing the outdated 
HYDROPLOT /HYDRLOG systems onboard its hydrographic vessels. 
The Shipboard Data System III (SDS III) is currently being 
built to increase automation and acquisition of hydrographic 
data. One of the capabilities of SDS III will be to acquire 
and process MLOP for sounding line positions. 

The methods presented in this thesis for data 
acquisition are based on the use of components, such as Del 
Norte Trisponders, which have been in existence since the 
early 1970' s. This thesis has demonstrated that using these 
systems in a MLOP configuration will substantially increase 
the accuracy of the survey. However, other systems exist 
which should be considered in conjunction with conventional 
MLOP. 

Krupp-Atlas Elektronik introduced a new short-range 
dynamic positioning system, Polarfix, in 1982. The 
advertised [Wentzell, 1983] advantages of the system are: 

(a) Range- Azimuth fixing to 0.1-m accuracy. 

(b) Automatic operator- free tracking of the vessel. 

(c) Portable shore stations. 
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Since the Polarfix system has the potential to alleviate 
some of the accuracy problems associated with the present 
two-LOP methods, additional studies should be made to 
ascertain the manufacturer’s claim as to the proport ed 
accuracy of the system. 

Differential GPS is a method which should be considered 
for large-scale hydrographic surveys. The theory has proven 
GPS extremely accurate (within centimeters) for control 
sites ashore where receivers acquire positional data from 
numerous passes of the satellites. Accuracies on the order 
of 2 to 3 m may be attained underway once the theory and 
technology have been developed. 

Future positioning systems may be integrated so that a 
combination of methods and systems can be analyzed to 
determine the ’’best" fix. 
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APPENDIX A 



FORTRAN PROGRAMS FOR DATA PROCESSING 

The following listed programs were used to process the 
LOP data for this thesis. Output files were generated to 
list the standard error of the adjusted coordinates, cy and 

cj, r ; the correlation coefficient, the semi-major and 

y 

semi-minor axes of the error ellipse, cj a and aj-, ; the 
orientation of the error ellipse, 0; and the standard error 
of position, a D • 

AZ2L0P 

Program AZ2L0P computes the positional accuracy of a fix 
obtained from two azimuth LOPs . This program is a 
modificaton of program UCOMPS . 

LSQAZ4 

Program LSQAZ4 computes positional accuracy using four 
azimuth LOPs. This program is a modification of program 
SILVA1 [Sliva, 1982] . The modification permits handling 
large amounts of data. 

RAZ2L0P 

Program RAZ2L0P computes positional accuracy of a fix 
obtained from one range LOP and one azimuth LOP. This 
program is a modificaton of program UCOMPS. It includes the 
error propagation algorithm. 
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RR+AZ 



Program RR+AZ computes the positional accuracy of a fix 
obtained from two range LOPs and one azimuth LOP. This 
program is a modification of program SILVA3 [Sliva, 1982] . 
The modification permits handling large amounts of data. 

RR2L0P 

Program RR2L0P computes positional accuracy using two 
range LOPs. This program is a modification of program 
UCOMPS . 

R2+AZ2 

Program R2+AZ2 computes the positional accuracy of a fix 
obtained from two range LOPs and two azimuth LOPs . This 
program is a modification of program SILVA3 [Silva, 1982]. 
The modification permits acceptance of another azimuth LOP. 

R3 + AZ 

Program R3+AZ computes the positional accuracy of a fix 
obtained from three range LOPs and one azimuth LOP. This 
program is a modification of program SILVA3 [Silva, 1982] 
which computed a least squares position from two range LOPs 
and one azimuth LOP. 

SLEAST 

Program SLEAST (Appendix B) computes the least squares 
position and positional accuracy of a fix from three or four 
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range LOPs. The program was written by Lt . Nick Perugini 
and the author from the algorithm presented by Cross [1981] . 

UCOMPS 

Program UCOMPS computes coordinates of a point in either 
Geographical Position (G.P.) or X and Y grid coordinates 
from range-range, range-azimuth, or azimuth- azimuth LOPs. 

It is a utility package program used by Hydrographic 
Sciences students at NPS. 
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FORTRAN PROGRAM FOR LEAST SQUARES POSITION FROM MULTIPLE RANGES 
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CORRECT OBSERVED RANGES FOR BASE-LINE CORRECTORS 
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FF=LANG-CCC 

DIFX=RANGE(1)*DC0S(FF) 
DIFY=RANGE ( 1 ) *DSIN ( FF ) 
FINX=XSTA( 1 ) -DIFX 



COMPUTE THE WEIGHT MATRIX 
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COMPUTE STANDARD DEVIATION OF X,Y, AND COVARIANCE = SIGX 
SIGY, AND SIGXY 

SIGX = SIGNO * DSQRT(Q(1, 1) ) 

SIGY = SIGNO * DSQRT (Q(2 , 2 ) ) 

SIGXY = ((SIGNO * SIGNO) * Q(l,2)) 
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PHI = (2.0D0 * Q( 1 , 2 ) ) / (Q(l,l) - Q(2,2)) 
THETA = 0.5D0 * DATAN (PHI) 

PI = DARCOS (-1.0D0) 

THETA = THETA * ( 180. 0D0 / PI) 
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PURPOSE: SOLUTION OF LINEAR LEAST SQUARES PROBLEM 

HIGH ACCURACY SOLUTION 
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9000 CONTINUE 

CALL UERTST (IER , 6HLLBQF 
9005 RETURN 
END 



APPENDIX C 



FOUR-RANGE MLOP FIX ACCURACIES 

EXAMPLE DATA OUTPUT FROM PROGRAM SLEAST 
POSITION FROM FOUR RANGES 



CONTROL 


STATIONS WERE 


SQUARE, 


CONK, 


USE MON 


, GEOCEIVER 


FIX# 


SIG X 


SIG Y 


SIG XY 


SIG A 


SIG B 


THETA 


SIGMA P 


1 


0,90 


1.13 


0.48 


1.43 


0.81 


-31.86 


1.45 


2 


1.34 


1.62 


0.89 


1.96 


1.22 


-32.67 


2. 10 


3 


1.39 


1.62 


0.76 


1.87 


1.28 


-32.87 


2.14 


4 


1.56 


1.79 


0.74 


1.98 


1.46 


-31.53 


2.37 


9 


0.31 


0.39 


0.00 


0.39 


0.31. 


-2.70 


0.50 


10 


0.76 


0.97 


-0.00 


0.97 


0.76 


0.30 


1.23 


11 


0.59 


0.76 


-0.01 


0.77 


0.59 


2.36 


0.96 


13 


0.71 


0.95 


-0.04 


0.96 


0.71 


5.69 


1.19 


14 


0.59 


0 , 80 


-0.04 


0.80 


0.59 


7.42 


0.99 


15 


0.59 


0,80 


-0.05 


0.81 


0.59 


9.04 


1.00 


17 


0.32 


0.45 


-0.02 


0.46 


0.32 


12.10 


0.56 


18 


0.47 


0.67 


-0.06 


0.69 


0.46 


13 . 39 


0.82 


20 


0.94 


1.39 


-0.31 


1.46 


0.92 


15 . 31 


1.68 


21 


1.09 


1.25 


-0.69 


1.67 


0.94 


37.36 


1.66 


22 


1.15 


1.28 


-0.64 


1.63 


1.01 


38.02 


1.72 


23 


1.70 


1.84 


-1.14 


2.23 


1.51 


38.71 


2.51 


24 


1.69 


1.79 


-0.92 


2.09 


1.52 


39.40 


2.46 
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FIX# 



SI G X SIG Y 



SIG XY SIG A 



SIG B THETA SIGMA P 



25 


0.98 


1.02 


- 0.25 


1.16 


0.89 


40.16 


1.41 


27 


1.12 


1.14 


- 0.20 


1.23 


1.04 


41.68 


1.59 


30 


0.89 


0.89 


- 0.05 


0.92 


0.86 


43.70 


1.26 


31 


0.93 


0.92 


- 0.03 


0.94 


0.91 


- 41.30 


1.31 


32 


0.19 


0.19 


- 0.00 


0.19 


0.19 


- 20.30 


0.27 


33 


0.86 


0.83 


0.01 


0.86 


0.83 


14.68 


1.20 


34 


0.79 


0.76 


0.02 


0.79 


0.75 


25.32 


1.09 


35 


0.96 


0.93 


0.05 


0.97 


0.92 


33.29 


1.34 


36 


0.88 


0.87 


0.05 


0.91 


0.85 


38.79 


1.24 


37 


0.87 


0.86 


0.06 


0.90 


0.83 


42.80 


1.22 


38 


0.57 


0.56 


0.03 


0.60 


0:54 


43.67 


0.80 


39 


0.87 


0.86 


0.10 


0.93 


0.82 


42.84 


1.23 


40 


0.50 


0.50 


0.04 


0.54 


0.46 


43.18 


0.70 


41 


0.58 


0.58 


0.06 


0.64 


0.54 


44 . 14 


0.82 


42 


0.72 


0.72 


0.10 


0.80 


0.65 


44.88 


1.01 


43 


0.67 


0.67 


0.10 


0.76 


0.60 


44.92 


0.94 


51 


0.32 


0.35 


0.02 


0.37 


0 ; 30 


- 32.03 


0.47 


53 


0.84 


0.93 


0.09 


0.96 


0.82 


- 22.93 


1.25 


54 


0.64 


0.72 


0.04 


0.73 


0.63 


- 17.40 


0.96 


56 


0.32 


0.37 


0.00 


0.37 


0.32 


- 5.92 


0.49 


57 


0.36 


0.42 


0.00 


0.42 


0.36 


- 2.39 


0.55 
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FIX# 



SI G X SIG Y 



SIG XY SIG A SIG B THETA SIGMA P 



58 


0.53 


0.63 


1 

o 

o 

o 


61 


0.51 


0.62 


- 0.02 


62 


0.45 


0.55 


- 0.02 


63 


0.51 


0.62 


- 0.03 


67 


1.80 


2.20 


- 0.71 


69 


1.06 


1.30 


- 0.31 


70 


0.29 


0.35 


- 0.03 


71 


0.36 


0.43 


- 0.05 


72 


0.81 


0.94 


- 0.27 


73 


1.14 


1.29 


- 0.63 


74 


0.83 


0.94 


- 0.39 


75 


0.75 


0.86 


- 0.39 



.63 


0.53 


1.44 


0.82 


.62 


0.51 


9.89 


0.80 


.56 


0.45 


11.76 


0.71 


.63 


0.51 


13.62 


0.81 


.30 


1.75 


20.89 


2.85 


.39 


1.02 


23.69 


1.67 


. 38 


0.28 


26.28 


0.46 


.48 


0.34 


30.20 


0.56 


.10 


0.74 


33.95 


1.24 


.62 


1.01 


36.58 


1.72 


.26 


0.71 


38.07 


1.25 


.28 


0.63 


38.51 


1.14 



0 

0 

0 

0 

2 

1 

0 

0 

1 

1 

1 

1 
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APPENDIX D 



FOUR-AZIMUTH MLOP FIX ACCURACIES 

DATA PROCESSED USING A MODIFICATION OF PROGRAM SILVA1 
FROM SILVA (1982) TO HANDLE LARGE DATA FILES. 

EXAMPLE DATA OUTPUT FROM PROGRAM LSQAZ4 . 

POSITION FROM FOUR AZIMUTHS. 

CONTROL STATIONS WERE SQUARE, CONK, USE MON, GEOCEIVER. 



FIX# 


SIG X 


SIG Y 


SIG XY 


SIG A 


SIG B 


THETA 


SIGMA 


76 


0.08 


0.15 


-0.00 


0.16 


0.08 


97.59 


0.17 


77 


0.04 


0.07 


0.00 


0.07 


0.04 . 


90.12 


0.08 


78 


0.61 


0.92 


0.05 


0.92 


0.61 


83.75- 


1.10 


79 


0.31 


0.41 


0.02 


0.41 


0.30 


78.10 


0.51 


80 ' 


0.23 


0.28 


0.01 


0.29 


0.23 


74.24 


0.37 


81 


0.80 


0.88 


0.05 


0.89 


0.79 


71.48 


1.19 


82 


0.43 


0.45 


0.00 


0.45 


0.43 


76.41 


0.62 


83 


0.42 


0.43 


-0.01 


0.43 


0.42 


129.29 


0.60 


84 


0.51 


0.51 


-0.02 


0.53 


0.49 


137.63 


0.72 


85 


0.25 


0.25 


-0.01 


0.27 


0.24 


135.85 


0.36 


86 


0.33 


0.33 


-0.02 


0.36 


0.31 


132.84 


0.47 


87 


0.39 


0.40 


-0.02 


0.42 


0.37 


131.27 


0.56 


88 


0.56 


0.55 


-0.03 


0.59 


0.53 


135.74 


0.78 


89 


0.43 


0.40 


-0.01 


0.44 


0.39 


158.24 


0.59 


90 


1.22 


0.99 


-0.09 


1.23 


0.99 


170.41 


1.57 
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FIX # 


SIG X 


SIG Y 


SIG XY 


SIG A 


SIG B 


THETA SIGMA P 


91 


0.58 


0.55 


- 0.10 


0.68 


0.50 


140.02 


0.80 


92 


0.20 


0.19 


- 0.01 


0.23 


0.18 


135.75 


0.28 


93 


0.73 


0.75 


- 0.12 


0.83 


0.67 


131.11 


1.04 


94 


0.58 


0.61 


- 0.07 


0.66 


0.54 


127.81 


0.84 


95 


0.61 


0.65 


- 0.07 


0.70 


0.58 


125.75 


0.89 


96 


0.44 


0.47 


- 0.03 


0.50 


0.42 


122.77 


0.65 


98 


0.33 


0.36 


- 0.01 


0.37 


0.32 


114.39 


0.49 


99 


0.51 


0.57 


- 0.02 


0.58 


0.50 


104.35 


0.76 


100 


0.39 


0.45 


- 0.01 


0.46 


0.39 


95.63 


0.60 


101 


0.39 


0.49 


- 0.00 


0.49 


0.39 


91.62 


0.63 


103 


0.76 


1.11 


- 0.04 


1.11 


• 0.76 


93.19 


1.34 


104 


0.20 


0.32 


- 0.01 


0.32 


0.20 


96.65 


0.38 


106 


0.32 


0.52 


- 0.06 


0.57 


0.31 


106.55 


0.61 


107 


0.28 


0.44 


- 0.03 


0.46 


0.27 


103.64 


0.52 


108 


0.70 


1.06 


- 0.13 


1.09 


0.69 


101.25 


1.27 


111 


0.61 


0.79 


- 0.05 


0.80 


0.60 


102.03 


0.99 


112 


0.54 


0.67 


- 0.05 


0.68 


0.53 


105.32 


0.86 


113 


1.47 


1.77 


- 0.39 


1.83 


1.44 


109.54 


2.30 


114 


0.53 


0.62 


- 0.06 


0.65 


0.51 


113.55 


0.81 


115 


1.00 


1.16 


- 0.24 


1.24 


0.96 


116.71 


1.54 


116 


0.60 


0.69 


- 0.10 


0.75 


0.57 


119.17 


0.91 


117 


0.74 


0.85 


- 0.17 


0.94 


0.69 


121.17 


1.13 


118 


0.71 


0.81 


- 0.17 


0.92 


0.66 


122.83 


1.08 


119 


0.56 


0.65 


- 0.12 


0.75 


0.51 


123.84 


0.86 


120 


0.45 


0.68 


- 0.12 


0.77 


0.42 


111.38 


0.81 
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FIX# SI G X 



SIG Y SIG XY SIG A 



SIG B THETA SIGMA P 



122 


1.44 


1.95 


- 0.97 


2.20 


1.37 


114.13 


2.43 


124 


0.41 


0.53 


- 0.06 


0.57 


0.39 


113.84 


0.66 


125 


0.53 


0.68 


- 0.10 


0.74 


0.51 


113.00 


0.86 


126 


0.67 


0.87 


- 0.14 


0.92 


0.65 


111.57 


1.10 


127 


0.50 


0.65 


- 0.07 


0.69 


0.49 


109.83 


0.82 


128 


0.51 


0.67 


- 0.07 


0.70 


0.50 


107.96 


0.85 


129 


0.45 


0.61 


- 0.05 


0.63 


0.44 


106.26 


0.75 


131 


0.36 


0.52 


- 0.04 


0.54 


0.35 


104.58 


0.63 


133 


1.03 


1.60 


- 0.47 


1.72 


1.00 


106.07 


1.90 


134 


0.31 


0.50 


- 0.05 


0.55 


0.30 


107 . 91 


0.58 


136 


0.44 


0.72 


- 0.14 


0.84 


0.42 


110.70 


0.85 


137 


0.61 


0.98 


- 0.23 


1.10 


0.59 


108.89 


1.15 


138 


0.76 


1.19 


- 0.30 


1.31 


0.74 


107.63 


1.42 


139 


0.79 


1.21 


- 0.28 


1.30 


0.77 


106.81 


1.44 


141 


0.27 


0.40 


- 0.03 


0.42 


0.26 


106.77 


0.48 


142 


0.34 


0.49 


- 0.04 


0.52 


0.33 


107.28 


0.60 


143 


0.51 


0.73 


- 0. 10 


0.77 


0.49 


107.89 


0.88 


145 


1.12 


1.63 


- 0.53 


1.76 


1.09 


108.65 


1.98 


147 


1.08 


1.67 


- 0.57 


1.82 


1.04 


107.68 


1.98 


148 


0.88 


1.45 


- 0.42 


1.58 


0.86 


106.30 


1.70 


150 


0.89 


1.30 


- 0.44 


1.48 


0.85 


112.26 


1.58 


151 


0.72 


0.92 


- 0.26 


1.09 


0.67 


119.01 


1.17 


152 


0.56 


0.62 


- 0.13 


0.75 


0.50 


127.65 


0.84 


154 


0.69 


0.60 


- 0.07 


0.73 


0.58 


155.12 


0.92 


155 


0.72 


0.62 


0.01 


0.72 


0.62 


3.15 


0.95 
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FIX # 


SIG X 


SIG Y 


SIG XY 


SIG A 


SIG B 


THETA ! 


SIGMA P 


156 


0.37 


0.26 


- 0.05 


0.47 


0.24 


151.08 


0.45 


157 


0.40 


0.38 


- 0.06 


0.51 


0.33 


138.38 


0.55 


158 


0.87 


0.92 


- 0.23 


1.06 


0.78 


129.72 


1.26 


159 


0.23 


0.25 


- 0.01 


0.26 


0.22 


124.69 


0.33 


160 


0.39 


0.40 


- 0.02 


0.43 


0.37 


127 . 75 


0.56 


161 


0.40 


0.41 


- 0.02 


0.44 


0.38 


131.07 


0.57 


163 


0.40 


0.42 


- 0.04 


0.46 


0.37 


127.67 


0.58 


164 


0.59 


0.65 


- 0.10 


0.72 


0.55 


123 . 37 


0.88 


165 


0.29 


0.35 


- 0.03 


0.38 


0.27 


118.55 


0.45 


166 


0.67 


0.86 


- 0.17 


0.95 


0.63 


114.46 


1.09 


167 


0.96 


1.70 


- 0.51 


1.82 


0.94 


103.91 


1.95 


168 


0.94 


1.55 


- 0.44 


1.67 


0.92 ' 


105.10 


1.81 


170 


1.64 


2.36 


- 1.05 


2.53 


1.60 


108.11 


2.88 


171 


0.78 


1.05 


- 0.20 


1.11 


0.76 


109.91 


1.30 


172 


1.01 


1.26 


- 0.28 


1.34 


0.97 


111.85 


1.62 


174 


1.81 


2.01 


- 0.57 


2.11 


1.75 


118.04 


2.71 


175 


0.23 


0.24 


- 0.01 


0.25 


0.22 


124.17 


0.33 


175 


0.23 


0.24 


- 0.01 


0.25 


0.22 


124.17 


0.33 


177 


1.07 


1.01 


- 0.11 


1.10 


0.98 


150.74 


1.47 


178 


1.41 


1.24 


- 0.17 


1.44 


1.23 


161.06 


1.88 


179 


0.22 


0.18 


- 0.00 


0.22 


0.18 


168.43 


0.28 
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PROGRAM " CHANGE3 FORTRAN Bl" WILL OUTPUT THE RANGES AND 
AZIMUTHS IN THE DESIRED ORDER AND ALSO OMIT BAD RANGES 
AND AZIMUTHS! ! ! ! 



POSITION FROM TWO RANGES AND TWO AZIMUTHS 
INTEGER N,I,J,K,IFIX,POS 

REAL* 8 PI , TB( 10 , 4 ) , E ,E1 , E2 , E3 , E4 ,XO , YO ,X01 , Y01 ,X02 , Y02 
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WRITE ( 3 , 150) 

150 FORMAT ( / / / , IX , 'FIX # ' , 2X , ' SIGO ' , 8X , ' SIG X',4X,'SI G Y’,4X, 
*' SIG P ' , 4X , ' SIG XY ' , 4X , ' SIG A’,3X,'SIG B 4X THETA ',/ ) 

IF SIGMA3=TB (3,4) IS IN METERS THE FOLLOWING IS COMMMENTED OUT 
IF SIGMA4=TB (4 , 4 ) IS IN METERS THE FOLLOWING IS COMMMENTED OUT 



CONVERSION DEGREES RADIANS (TB (3 , 4 ) ) 
TB(3,4)=TB(3,4)*(PI/180.0) 
TB(4,4)=TB(3,4)*(PI/ 180.0) 
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GO TO 60 
CONTINUE 

OMEGA=DATAN(2.0D0*Q(l,2)/ (Q(l , 1)-Q(2 ,2 ) ) ) 
IF (OMEGA. LT.O.ODO)GO TO 61 
GAMA=OMEGA/ 20 . 0D0 
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APPENDIX F 



TWO- RANGE AND TWO -AZIMUTH MLOP FIX ACCURACIES 
EXAMPLE DATA OUTPUT FROM PROGRAM R2+AZ2 
POSITION FROM 2 RANGES AND 2 AZIMUTHS 

CONTROL STATIONS WERE SQUARE, SQUARE ECC . , USE MON, USE MON 



FIX# 


SIG X 


SIG Y 


SIG XY 


SIG A 


SIG B 


THETA ; 


SIGMA P 


1 


1.47 


1.66 


0.28 


1.71 


1.44 


68 . 63 


2.22 


2 


1.81 


2.04 


0.42 


2.09 


1.78 


67.92 


2.73 


4 


1.97 


2.16 


0.45 


2.22 


1.93 


65.66 


2.93 


6 


0.97 


1.04 


0.10 


1.07 


0.95 


63.17 


1.43 


7 


0.84 


0.89 


0.06 


0.91 


0.82 


61.98 


1.22 


8 


0.57 


0.60 


0.02 


0.61 


0.56 


60.74 


0.82 


9 


0.85 


0.88 


0.05 


0.90 


0.84 


59.44 


1.23 


10 


0.99 


1.01 


0.05 


1.03 


0.97 


58 . 19 


1.41 


11 


0.71 


0.72 


0.02 


0.73 


0.70 


56.92 


1.02 


13 


0.71 


0.71 


0.00 


0.71 


0.71 


147.31 


1.01 


14 


0.40 


0.40 


-0.00 


0.41 


0.40 


143.00 


0.57 


15 


0.83 


0.82 


-0.03 


0.85 


0.81 


141.40 


1.17 


17 


0.58 


0.57 


-0.03 


0.60 


0.55 


138.59 


0.81 


18 


0.55 


0.55 


-0.04 


0.58 


0.52 


137.41 


0.78 


20 


0.24 


0.24 


-0.01 


0.26 


0.22 


135.06 


0.33 


21 


1.00 


1.11 


-0.18 


1.17 


0.96 


118.66 


1.49 


22 


0.87 


0.96 


-0.13 


1.01 


0.84 


118.64 


1.30 
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FIX # 


SIG X 


SIG Y 


SIG XY 


SIG A 


SIG B 


THETA SIGMA P 


23 


1.43 


1.55 


- 0.29 


1.61 


1.38 


118.71 


2.11 


24 


1.22 


1.31 


- 0.18 


1.36 


1.19 


118.96 


1.79 


25 


1.04 


1.10 


- 0.11 


1.13 


1.01 


119.51 


1.51 


27 


0.53 


0.55 


- 0.02 


0.56 


0.52 


121.31 


0.76 


30 


1.10 


1.11 


- 0.04 


1.13 


1.09 


125.16 


1.56 


31 


0.64 


0.65 


- 0.01 


0.65 


0.64 


126.41 


0.91 


32 


0.72 


0.72 


- 0.00 


0.72 


0.72 


126.73 


1.02 


33 


0.79 


0.78 


0.01 


0.79 


0.78 


38.45 


1.11 


34 


0.85 


0.84 


0.02 


0.86 


0.83 


39.87 


1.19 


35 


0.64 


0.64 


0.01 


0.65 


0.63 


41.54 


0.91 


.36 


1.11 


1.11 


0.05 


1.14 


1.09 


43.24 


1.57 


37 


0.52 


0.52 


0.01 


0.53 


0.51 


44.92 


0.74 


38 


0.69 


0.69 


0.03 


0.71 


0.68 


46 . 49 


0.98 


39 


0.76 


0.77 


0.04 


0.80 


0.74 


47.84 


1.09 


40 


0.72 


0.73 


0.04 


0.76 


0.70 


49 . 30 


1.03 


41 


0.65 


0.66 


0.04 


0.69 


0.63 


50.80 


0.93 


42 


1.04 


1.06 


0.11 


1.10 


1.00 


52.39 


1.48 


43 


0.99 


1.02 


0.10 


1.06 


0.95 


53.78 


1.42 


44 


1.38 


1.44 


0.22 


1.50 


1.33 


55.09 


1.99 


47 


1.69 


1.84 


0.38 


1.91 


1.65 


61.79 


2.50 


48 


2.34 


2.53 


0.70 


2.62 


2.28 


61.34 


3.45 


49 


2.08 


2.23 


0.52 


2.31 


2.02 


61.06 


3.05 


51 


0.97 


1.02 


0.10 


1.05 


0.94 


59 o 55 


1.41 


52 


1.16 


1.21 


0.13 


1.25 


1.13 


58.34 


1.68 


53 


1.41 


1.47 


0.17 


1.51 


1.38 


57.07 


2.03 
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54 

56 

57 

58 

61 

62 

63 

67 

69 

70 

71 

73 

74 

75 



SIG X SIG Y SIG XY SIG A SIG B THETA SIGMA P 



1.32 


1.35 


0.13 


1.39 


1.29 


55.64 


1.89 


0.72 


0.73 


0.02 


0.75 


0.71 


52.97 


1.03 


0.74 


0.74 


0.02 


0.75 


0.73 


51.29 


1.05 


0.87 


0.87 


0.01 


0.87 


0.86 


50.01 


1.23 


1.18 


1.18 


- 0.06 


1.21 


1.16 


136 . 12 


1.67 


0.28 


0.28 


- 0.00 


0.29 


0.28 


134.85 


0.40 


0.76 


0.76 


- 0.04 


0.79 


0.74 


133.41 


1.08 


1.14 


1.18 


- 0.18 


1.25 


1.09 


128.73 


1.64 


0.82 


0.86 


- 0.12 


0.92 


0.77 


127.53 


1.19 


0.32 


0.34 


- 0.02 


0.37 


0.31 


126.80 


0.47 


0.45 


0.49 


- 0.04 


0.53 


0.43 


125.90 


0.67 


0.83 


0.90 


- 0.16 


0.99 


0.78 


124.70 


1.23 


0.94 


1.03 


- 0.22 


1.13 


0.88 


124.25 


1.40 


0.92 


1.00 


- 0.22 


1.11 


0.85 


124.50 


1.36 



118 



LIST OF REFERENCES 



Anonymous, Microwave Trisponder System Handbook , Document No. 
1400072C, Del Norte Technology, Inc., 1982. 

Anonymous, International Hydrographic Organization Special 
Publication No. 44 (2nd Edition), International Hydrographic 
Organization Standards for Hydrographic Surveys and Classification 
Criteria for Deep Sea Soundings , 1982. 

Bomford, G., Geodesy , Fourth Edition, Clarendon Press, 1980. 

Casey, M. J., Multiranging - An Expensive Solution to a Simple 
Problem , Lighthouse, Ed. No. 25, April 1982. 

Cross, P. A., The Computation of Position at Sea , The 
Hydrographic Journal, No. 20, April 1981. 

Davis, R. E. , and others, Surveying : Theory and Practice , 

Sixth Edition, McGraw-Hill, 1981. 

Eaton, R. M. , Smoothing and Multi-Ranging , Lighthouse, Ed. No. 
25, April 1982. 

Greenwalt , C. R. , and M. Schultz, AQIC Technical Report No. 
96, USAF Aeronautical Chart and information Center, St. 
Louis, Mo., Principals of Error Theory and - Cartographic 
Applications , 1962. 

Heinzen, M. R. , Hydrographic Surveys : Geodetic Control Criteria , 

M.S. Thesis, Cornell University, 1977. 

Holder, R. H. , National Ocean Service, Atlantic Marine 
Center, Operations Order 79, Del Norte Short Range Electronic 
Control System, 1983. 

Ingham, A. E., Hydrography for the Surveyor and Engineer, Second 
Edition, Granada Publishing Ltd., 1984. 

Kaplan, A. , Error Analysis of Hydrographic Positioning and the 
Application of Least Squares, M.S. Thesis, Naval Postgraduate 
School, Monterey, 1980. 

Mikhail, E. M. , Observations and Least Squares, IEP , 1976. 

Mikhail, E. M. , and G. Grade, Analysis and Adjustment of 
Survey Measurements, Von Nostrand Reinhold, 1981. 



119 



Mueller, I. I., and K. Ramsayer, Introduction to Surveying , 
Frederick Unger Publishing Co., 1979. 

Munson, R. C. , Positioning Systems : Report on the Work of WG4l4b , 

paper presented at the XV International Congress of 
Surveyors, Stockholm, Sweden, 1977. 

Perugini , N. E. , Criteria for the Classification of Hydrographic 
Positioning Data , M.S. Thesis, Naval Postgraduate School, 
Monterey, 1984. 

Riemersma, Ir. J. G. , Quality Control of Offshore Positioning 
Systems , International Hydrographic Review, LVII (2), July 
1980. 

Saxena, N. K. , Improvement of a Geodetic Triangulation Through 
Control Points Established by Means of Satellites or Precision 
Traversing , Report No. 177, The Ohio State University 
Research Foundation, 1972. 

Saxena, N. K. , Applications of Marine Geodesy in Support of 
National Objectives in Ocean Science , Engineering and Operations , 
Marine Technology Society, 1980. 

Silva, F. C. , Calculation of Hydrographic Position Data by Least 
Squares Adjustment , M.S. Thesis, Naval Postgraduate School, 
Monterey, 1982. 

r 

Umbach, M. J., Hydrographic Manual Fourth Edition , National 
Oceanic and Atmospheric Administration, National Ocean 
Survey, 1976. 

Wallace, J. L. , National Ocean Survey Technical Manual 
No. 1., HYDROPLOT I HYDROLOG System Manual , 1971. 

Waltz, D. A. , An Error Analysis of Range- Azimuth Positioning , 

M.S. Thesis, Naval Postgraduate School, Monterey, 1983. 

Weeks, C. , Some thoughts on the processing of multiple posit ion 
lines and on the vouge for least squares adjustment , The 
Hydrographic Journal, No. 23, January 1982. 

Wentzell, H. F. , Range / Azimuth Position Fixing and Data Processing 
System with Decimeter Accuracy for Civil Engineer ing Survey , The 
Hydrographic Journal, No. 27, January 1983. 

Whitsell, G. R. , and K. Berstis, Del Norte DDMU 520 Short Range 
Positioning System Stability and Range Accuracy Test Results , 

Systems Planning and Applied Technology Staff, National 
Ocean Service, 1983. 



120 



INITIAL DISTRIBUTION LIST 



No. 



1. Defense Technical Information Center 
Cameron Station 

Alexandria, VA 22304-6143 

2. Library, Code 0142 
Naval Postgraduate School 
Monterey, CA . 93943-5002 

3. Chairman, Department of Oceanography 
Code 68Mr 

Naval Postgraduate School 
Monterey, CA 93943-5000 

4. N0AA Liaison Officer 
Post Office Box 8688 
Monterey, CA 93943-0688 

5. Professor Narendra Saxena 
Department of Civil Engineering 
University of Hawaii at Manoa 
336 Holmes Hall, 2540 Doles St. 

Honolulu, Hawaii 96822 

6. Lt . Samuel P. De Bow, Jr. 

N/M01, Room 214, WSC-1 

National Oceanic and Atmospheric Administration 
Rockville, MD 20852 

7. Office of the Director 

Naval Oceanography Division (OP-952) 

Department of the Navy 
Washington, D.C. 20350 

8 . Commander 

Naval Oceanography Command 
NSTL , MS 39529 

9. Commanding Officer 

Naval Oceanographic Office 
Bay St . Louis 
NSTL, MS 39522 



Copies 

2 

2 

1 

1 

1 

1 

1 

1 

1 



121 



10 . 



1 



Commanding Officer 

Naval Ocean Research and Development Activity 
Bay St. Louis 
NSTL , MS 39522 

11. Chief of Naval Research 1 

800 N. Quincy Street 

Arlington, VA 22217 

12. Chairman, Oceanography Department 1 

U. S. Naval Academy 

Annapolis, MD 21402 

13. Chief, Hydrographic Programs Division 1 

Defense Mapping Agency (Code PPH) 

Building 56 

U.S. Naval Observatory 
Washington, D.C. 20305-3000 

14. Associate Deputy Director for Hydrography 1 

Defense Mapping Agency (Code DH) 

Building 56 

U.S. Naval Observatory 
Washington, D.C. 20305-3000 

15. Topographic Sciences Department 1 

Code TSD-MC 

Defense Mapping School 
Ft. Belvoir , VA 22060-5828 

16. Director, Charting and Geodetic Services 1 

N/CG, Room 1006, WSC-1 

National Oceanic and Atmospheric Administration 
Rockville, MD 20852 

17. Chief, Nautical Charting Division 1 

N/CG2, Room 1026, WSC-1 

National Oceanic and Atmospheric Adminsitrat ion 
Rockville, MD 20852 

18. Chief, Hydrographic Surveys Branch 1 

N/CG24, Room 404, WSC-1 

National Oceanic and Atmospheric Administration 
Rockville, MD 20852 

19. Program Planning, Liaison, and Training Div. 1 

NC2, Room 105, Rockwall Building 

National Oceanic and Atmospheric Administration 
Rockville, MD 20852 



122 



20 . 



1 



Director, Pacific Marine Center 
N/MOP 

National Ocean Service, NOAA 
1801 Fairview Avenue, East 
Seattle j WA 98102 

21. Director, Atlantic Marine Center 1 

N/MOA 

National Ocean Service, NOAA 
439 West York Street 
Norfolk, VA 23510 

22. Chief, Program Services Division 1 

N/MOA2 

National Ocean Service, NOAA 
439 West York Street 
Norfolk, VA 23510 

23. Mr. Samuel P. De Bow, Sr. 1 

2622 Ann Street 

Philadelphia, PA 19134 

24. Lt . Nicholas E. Perugini 1 

NOAA Ship RAINIER 

1801 Fairview Avenue, East 
Seattle, WA 98102 

25. Mr. Jon Didrick Andersen 1 

Norges Sjokartverk 

Post Boks 60 
N-4001 Stavanger 
NORWAY 

26. Mr. Kenneth G. Vadnais 1 

6252 Covington Way 

Goleta , CA 93117 

27. Mr. Albert Visek 1 

1306 Hartranft Rd. 

Fort Washington, PA 19034 

28. CDR Thomas W. Ruszala 1 

NOAA Officer Training Center 

c/o U.S. Merchant Marine Academy 
Furuseth Hall, 3rd Deck 
Kings Point, NY 11024 



123 






Thesis 
D18884 
c. 1 



/ 




De Bow 

Application of multi 
pie lines of position 
for hydrography. 




) * r j O ", 

'„IOv « 

Thesis 

D13834 Bow 

Application of multi 

pie lines of position 

for hydrography. 



thesD 18884 

Application of multiple lines of positio 




3 2768 000 65813 2 

DUDLEY KNOX LIBRARY 






